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Abstract 

We calculate the transverse-momentum Qt spectrum of the dilepton in the trans- 
versely polarized Drell-Yan process on the basis of the factorization theorem in QCD. 
We take into account universal logarithmically enhanced corrections in edge region of 
phase space by resumming multiple soft-gluon emissions to all orders in the small Qt 
region. 



*) Deceased. 
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§1. Introduction 



The nucleoli structure appearing in high energy processes has been studied for a long time 
on the basis of quantum chromodynamics (QCD). The basic theoretical framework has been 
developed as the "factorization theorem" in QCD. 1 ) As a result of factorization, the physical 
quantities (cross sections) are obtained as a convolution of the short- and long-distance 
parts. The former contains all the dependence on the hard scale, and the latter is controlled 
by the nonperturbative dynamics of QCD. We can apply perturbation theory to the short- 
distance part, thanks to the asymptotically free nature of QCD. The long-distance parts, 
however, can be determined only by experiments or some nonperturbative method, such as 
lattice QCD. The importance and advantage of perturbative QCD based on the factorization 
theorem reside in the fact that this theorem allows us to define the long-distance parts 
as process-independent universal objects, which are represented unambiguously as nucleon 
matrix elements of the operators of quarks and/or gluons and are often much simpler than 
the original quantities. 

For spin-independent processes, the many experiments performed to this time have helped 
us to determine the long-distance, nonperturbative parts as "parton distributions" inside the 
nucleon, 2 ^ and we have obtained a consistent understanding of the perturbative and nonper- 
turbative dynamics, except in "edge regions" of phase space. ^ For spin-dependent processes, 
however, our understanding is still poor, and many questions remain unanswered. 3 )' 4 ) There- 
fore, to understand spin-dependent processes and the spin structure of nucleons through them 
is an important problem. Furthermore, spin-dependent quantities are, in general, believed 
to be quite sensitive to the structure of interactions among particles. These are the reasons 
for the great amount of activity recently in high-energy spin physics. 

It is expected that a number of ongoing polarization experiments, such as RHIC-Spin, 5 )' 6 ) 
using polarized proton-proton collisions, 7 ) HERMES 8 ) and COMPASS, 9 ) using lepton scat- 
tering off polarized protons, etc., will provide important experimental data to reveal spin- 
dependent phenomena associated with the structure of nucleons. We also expect that data 
from future polarization experiments using proton-proton collisions at J-PARC, 10 ) proton- 
antiproton collisions at GSI, 11 ) etc will provide useful information. Therefore, it is important 
and interesting to investigate various processes to be studied in those experiments. 

Among the many spin-dependent processes, the polarized Drell-Yan (DY) process plays 
a unique role. Indeed, transversity 12 ) -14 ) is one of the characteristic observables which can 
be measured in DY using a transversely polarized beam. The transversity 5q(x) is a twist- 
2 parton distribution associated with the probability distribution of transversely polarized 
quarks inside transversely polarized nucleons, i.e., the partonic structure of nucleons which 
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is complementary to that represented by the other twist-2 distributions, such as the familiar 
density and helicity distributions q(x) and Aq(x). However, 8q(x) is not yet well understood, 
because 5q(x) cannot be measured in inclusive deep inelastic scattering (DIS) due to its 
"chiral-odd" nature. 12 )~ 14 ) In this paper, we focus on the transversely polarized Drell-Yan 
(tDY) process, p^p^ — > l + l~X, producing the dilepton l + l~ with an invariant mass Q ^> /Iqcd 
in the final state, as a process that is likely to allow investigation of the transversity Sq(x). 
Our aim is to calculate the QCD corrections to the tDY cross section on the basis of the 
QCD factorization framework, establishing control over the large higher-order corrections 
near edge region of phase space. 

Simpleminded calculations of QCD corrections to the DY cross section suffer from ultra- 
violet (UV) and infrared (IR) divergences due to the loops associated with (massless) quarks 
and gluons. The UV divergence can be regularized and renormalized straightforwardly using 
a standard procedure, and therefore it does not pose any problem. By contrast, the IR diver- 
gence is more intricate and must be treated using the factorization theorem: x ) Introducing 
an appropriate IR regulator, one has to confirm that the IR divergences from the different di- 
agrams cancel ( "Kinoshita-Lee-Nauenberg cancellation" 15 )) or are completely factorized into 
parton distribution functions. Such a calculation has been performed for 0(a s ) corrections 
to tDY using various IR regularization schemes. For example, the one-loop calculation of 
the tDY cross section was done by Vogelsang and Weber 16 ) using the massive gluon scheme. 
The same calculation was also done using the dimensional reduction scheme. 17 ' The relations 
among the results obtained with different schemes are discussed in Ref. 18). The result in 
the dimensional regularization scheme was obtained 19 -* by using the scheme transformation 
relation (see also Ref. 20)). We also note that all these works 16 - 1-20 ) investigating tDY treat 
the case in which the transverse momentum Qt of the final dilepton with respect to 
the beam axis is unobserved (integrated).*) 

When a calculation of such QCD corrections is performed, one encounters an interesting 
technical problem in the case of transverse polarization. As is well known, 12 ) the cross 
section depends on the azimuthal angle, 0, of the observed particle, and its dependence in 
tDY is of the form cos(20). Therefore, we must keep the azimuthal angle dependence of the 
cross section in the case of transverse polarization, which makes it difficult to perform the 
phase space integrals in the higher-order calculations. This difficulty becomes much more 
severe if dimensional regularization is employed. One cannot use the techniques developed 
for the unpolarized and longitudinally polarized DY. A related problem in the dimensional 
regularization scheme has been discussed by Kamal, 21 ) but an explicit result was not given. 

*) Ref. 16) treats also the case with Qt observed, restricting Qt to large values, which makes the IR 
regulator irrelevant for calculation of the 0(a s ) contributions to the cross section. 
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Recently, Mukherjee et al. ' proposed a new technique which allows one to overcome a 
similar problem in their calculation of prompt photon production. 

In this paper, we employ the dimensional regularization scheme and discuss our approach 
to directly integrate out the phase space in D dimensions. The final result for the tDY cross 
section up to 0(a s ) accuracy in the MS scheme has been reported in a previous paper, 23 -* and 
here we present the details of our calculation. We believe that our explicit calculation in the 
dimensional regularization is interesting and useful. Once the above-mentioned complication 
associated with the transverse polarization is worked out, the dimensional regularization pro- 
vides us with a manifestly gauge-invariant and the most transparent framework to calculate 
the QCD corrections for both the Qr-unobserved and Qr-observed cases. For the latter case, 
in particular, we are able to derive the tDY cross section up to 0(a s ) without any restriction 
on Qt, explicitly isolating the terms that diverge as ln(Q 2 /Qt)/Qti ^-/Qti anc ^ &{Qt)i as 
Qt — > 0. The resulting "QT-differential" cross section gives the leading order (LO) QCD 
prediction in the large Qt region, in which Qt ~ Q, where the lepton-pair production via 
the DY mechanism has to be accompanied by the radiation of at least one recoiling "hard" 
gluon, and for this reason the fixed-order truncation of perturbation theory is effective. On 
the other hand, the unlimited growth of the singular terms of the form ln(Q 2 /Q T )/Q T in the 
cross section at small Qt is associated with the recoil from "soft" gluon radiation, and this 
implies that we cannot truncate perturbation theory and have to calculate the higher-order 
QCD corrections beyond the one-loop level in this edge region of phase space. 

The higher-order corrections that control the small Qt behavior of the cross section 
can be taken into account through the soft gluon resummation ("transverse-momentum 
(Qt) resummation"). As the next step, we derive the soft gluon resummation for our Qt- 
differential tDY cross section, so that we can extend our LO QCD prediction at large Qt to 
the entire range of Qt- Below, we summarize the development of the soft gluon resummation. 

A small value of Qt (<C Q) implies that there exists a new scale in the problem, and as 
a result, in the perturbative calculation, there appear terms containing large logarithms of 
Q 2 / Qt'- The coefficient of a™(Q 2 ) includes a factor of 1/ Q T multiplying a series of logarithms 
of the form \n m (Q 2 /Q T ), with m = 0, 1, . . . , 2n — 1. The pattern of these logarithmic 
terms is characteristic of a theory with massless vector bosons, such as QCD and QED, 
and is produced by recoil from the radiation of gluons and photons. The first work dealing 
with these enhanced "recoil logarithms" in QCD was carried out by Dokshitzer, Dyakonov 
and Troyan. 24 ) Their result corresponds to the leading logarithmic (LL) resummation in 
momentum space, i.e., the resummation of the terms a™ In 211-1 (Q 2 / Q 2 ^) to all orders in a s . 
The level of this approximation is the same as that in Sudakov's QED analysis, 25 ^ and the 
result was derived by imposing the "strong ordering" of the gluon's transverse momenta 
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in the relevant Feynman diagrams; the strong ordering excessively constrains the phase 
space of the emitted soft gluons, and thus results in the transverse-momentum conservation 
being ignored. To take into account the transverse-momentum conservation, Parisi and 
Petronzio 26 ** developed a formulation in the space of the impact parameter, b, which is 
the Fourier conjugate of the Qt space. The relation between the Qr-space and the b- 
space approaches is analyzed in Ref. 27). That work clarifies the impact of the transverse- 
momentum conservation on the subleading logarithmic terms. The general form of the 
soft gluon resummation to all orders of logarithms was proved and formulated in the b 
space approach by Collins and Soper. 28 ) The universal two-loop anomalous dimension in 
the resummed cross sections, which is necessary for the next-to-leading logarithmic (NLL) 
analysis of all relevant processes, was first calculated in Ref. 29). This result was based on a 
plausible assumption, and Davies et al. 30 -* confirmed the result using an explicit calculation 
of the unpolarized DY process up to order a 2 . They also carried out a phenomeno logical 
analysis, but only for small Qt- 

Advanced 6-space formulations of the soft gluon resummation, which are suitable for 
phenomenological analyses of DY in all Qt regions, were developed in the following works. 
Altarelli et al. 31 -* proposed a recipe to include the NLL resummation effects into the un- 
polarized DY cross section, but their approach is somewhat naive; an application of this 
approach to the longitudinally polarized DY process is considered in Ref. 33). Presently, 
the formulation which is valid to all orders of logarithms developed by Collins, Soper and 
Sterman (CSS) 32 ) is regarded as standard. Recently, de Florian and Grazzini 34 ) derived a 
universal expression for the CSS resummation formula up to next-to-next-to-leading logarith- 
mic (NNLL) accuracy, which is applicable to DY production, electroweak boson production, 
Higgs boson production, etc. Further developments of formulations and their applications 
to phenomenology are underway (see Refs. 35)-37) and references therein). *) The extension 
of the formulations to DY with transversely polarized beams has been carried out by the 
present authors, and the main results are reported in Refs. 23) and 56).**- 1 

The general formulation of soft gluon resummation used to obtain our results is only 
briefly described in Refs. 23) and 56). In this paper, we discuss the general form of the 
6-space resummation formula in detail from a modern viewpoint, emphasizing its theoretical 
basis as well as its physical content. We also demonstrate that the b space can be divided 

*) Besides the Qt resummation discussed in this paper, another kind of soft gluon resummation, the 
so-called "threshold resummation," has been developed for the purpose of resumming the large higher-order 
corrections near the threshold of partonic scattering. 41 )' 42 )' 53 )' 68 ) Also, the joint resummation formalism 
has been devised to perform the Qt and threshold resummations simultaneously. 45 )~ 47 ) 

**) Sec Ref. 38) for a treatment within the LL-level resummation. The extension to polarized semi- 
inclusive DIS has recently been discussed. 39 ' 
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into three distinct regions, associated with different distance scales; each of these three 
regions has to be treated differently, and the results for the three regions must eventually 
be combined in a consistent manner. This point is addressed in the original papers on the 
6-space approach, 26 )' 32 ) but here we present a more systematic explanation of this point, 
which is important to ensure the maximal applicability of the resummation formalisms to a 
wide range of processes, including spin-dependent processes at moderate and high energies. 
Then, we apply the resummation formalism to our tDY QT-differential cross section, which 
allows us to include all orders resummation of the logarithmically-enhanced contributions 
for small Qt due to multiple emissions of soft gluons in QCD. We perform the corresponding 
resummation up to NLL accuracy, and the result is combined with the fixed-order LO cross 
section that controls the large Qt region, yielding a tDY cross section with uniform accuracy 
over the entire range of Qt- We also explicitly derive the important properties of our cross 
section, which were briefly addressed in previous papers. 23 -*' 56 - ) For example, our cross section 
satisfies the unitarity constraint exactly; the NLL resummation is controlled completely by 
the "saddle point" in the b space in the asymptotic regime, Q ^> Aqcb, Qt ~ 0. As 
an application of our results, we calculate the dilepton Q T spectrum and the cross section 
asymmetry as functions of Qt- These quantities are to be observed in polarized pp collisions 
with large CM energy > 200 GeV at RHIC. We also present the results to be observed 
in polarized pp collisions with moderate \/S = 10 GeV at J-PARC. We demonstrate that, 
for both RHIC and J-PARC energies, the soft gluon resummation is crucial for making a 
reliable QCD prediction for the small Qt region, where the bulk of dileptons is produced. 

To our great sorrow, one of the present authors, Jiro Kodaira, died on September 16, 
2006. The present work was performed by the three authors jointly, and many parts of this 
paper are based on our "collaboration notes," originally written by Jiro. To complete this 
paper, we have had to reorganize and expand those notes without the direct assistance of 
Jiro, but we have been guided by his style in the approach to the problem, which we learned 
through our association with him for more than ten years. 

The rest of the paper is organized as follows. Sections 2-4 are devoted to the calculation 
of the tDY cross section up to 0(a s ) in the dimensional regularization. Section 2 is mainly 
introductory, explaining the factorization theorem for the tDY cross section, and the relevant 
partonic mechanism to 0(a s ). The total and QT-differential cross sections are derived in §§3 
and 4, respectively, by performing a collinear factorization in the MS scheme. Sections 5-9 
contain discussion of the soft gluon resummation. In §5 we introduce the general formalism of 
soft gluon resummation, and in §6 we elaborate on its 6-space structure. We perform the soft 
gluon resummation for tDY up to NLL accuracy in §7. In §8, we also derive an asymptotic 
formula for our NLL resummed cross section in the Q T ~ region. Section 9 contains 
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numerical results for tDY processes to be observed at RHIC and J-PARC. We demonstrate 
the roles of QCD soft gluon effects in the cross sections and asymmetries as functions of 
Qt, and discuss a new approach to extracting the transversity from experimental data. 
The final section, §10, is reserved for conclusions. This paper contains three appendices. 
In Appendix A, we collect the operator definitions and basic properties of transversity, in 
Appendix B we collect the formulae for the tDY cross sections integrated over the rapidity 
of the dilepton, and Appendix C supplements the discussion given in §7. 

§2. Drell-Yan mechanism to order ot s 

2.1. Factorization formula for the transversely polarized Drell-Yan process 
The process we consider is the tDY process, 

h ± (Pi , Si) + h 2 {P 2 , S 2 ) - 1(h) + l(k 2 ) + X, (2-1) 

where hi and h 2 denote spin-1/2 hadrons with momenta Pi and P 2 and transverse spins Si 
and S 2 , and Q = k 1 + k 2 is the 4-momentum of the DY pair. The tDY (2-1) may be induced 
by partonic subprocesses, such as quark-antiquark annihilation, 

q(Pu si) + q(p2, s 2 ) -> 1(h) + l(k 2 ) + X, (2-2) 

where q and q have momenta pi and p 2 and transverse spins si and s 2 . This can be formulated 
precisely using QCD factorization, 1 ^ which is valid when Q 2 is large, i.e., Q 2 3> Aq CD , w ^ n 
r = Q 2 /S (S = (Pi+P 2 ) 2 ) fixed. In this case, (2-2) and q( Pl , Si)+q(p 2 , s 2 ) -> l(ki)+l(k 2 )+X 
provide the dominant mechanism for tDY in the region with pi = X{P{ (i = 1,2), where x\ and 
x 2 denote the longitudinal-momentum fractions of partons inside the parent hadrons hi and 
h 2 , respectively, while the contributions from the other regions or other partonic processes 
give subleading corrections. The values that xi or x 2 actually takes are determined by the 
parton distributions, due to the long-distance, nonperturbative dynamics inside each hadron 
/ii 5 2, whose typical distance scale ~ 1/Mqcd i s well-separated from the short-distance scale 
of order 1/Q, within which partonic annihilation processes like (2-2) occur. Therefore, the 
short- and long-distance contributions can be treated separately, and the spin-dependent 
cross section, A T da = (da (Si, S 2 ) — da(Si, —S 2 ))/2, for tDY (2-1) is given as a convolution 
of the short- and long-distance contributions, over each of the partonic momentum fractions 
Xi and x 2 for the two hadrons hi t2 , 12 ^ 13 ^ 

A T da = / dx 1 dx 2 5H(xi, x 2 ;fi 2 F ) [A T da(p 1 ,p 2 ;Q 2 ; /40] Pi=!B .p. > ( 2 ' 3 ) 
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up to the higher-twist corrections suppressed by the powers of 1/Q. Here, fi F is the fac- 
torization scale, AtcI(j(pi,P2; Q 2 ', h 2 f ) is the transversely polarized partonic cross section for 
(2-2) representing the short-distance contribution, and 

SH(x 1 ,x 2 ;/j 2 f ) = ^24 [Sq hl (x 1: ii F )Sq h2 (x2, H 2 F ) + Sq hl (x 1: /i 2 F )Sq h2 (x 2: ii 2 F )] (2-4) 
i 

is the long-distance contribution as a product of the transversity distributions 12 )' Sq hl (x 1: jip) 
and 5qh 2 (x 2 , [ip) for the two colliding hadrons, hi and h 2 , summed over the massless quark 
flavors q with their charge squared e 2 . 

Intuitively, the leading-twist factorization formula (2-3) can be interpreted using the lan- 
guage of the parton model, similarly to the factorization formula for the twist- 2 structure 
functions in DIS. In particular, the transversity distribution 5qh(x,/i F ) is one of the three 
independent twist- 2 quark distribution functions for a spin- 1/2 hadron h, which are asso- 
ciated with the probability distribution of finding a quark with a longitudinal momentum 
fraction x inside the parent hadron h. In fact, 5qh(x,fi F ) is defined by 12 ^ 14 ^ 

5q h (x, /4) = q[(x, /4) - q l h (x, /4) , 

where q\ (qjj denotes the probability distribution of transversely polarized quarks with spin 
parallel (antiparallel) to the spin of their transversely-polarized parent hadron. 14 ) The spin- 
dependent partonic cross section participating in (2-3) is defined as 

X 

A T da(p 1 ,p 2 ;Q 2 ;fi 2 F ) = - [da(p 1 ,p 2 ; si , s 2 ;Q 2 ; fi 2 F ) 

- da(p 1 ,p 2 ;s 1 , -s 2 ;Q 2 ;n 2 F )] , (2-5) 

in terms of the partonic hard cross section da(pi,p 2 ; Si , s 2 ;Q 2 ;fi 2 F ) which describes the 
annihilation process (2-2) of the quark and antiquark with the transverse spins Si and s 2 , 
respectively. Note that at the leading twist level, the gluon distributions do not contribute to 
(2-3) for the transversely polarized process, due to its chiral-odd nature. 12 ^ 14 -* Therefore, only 
the qq annihilation subprocesses give the relevant twist- 2 mechanism for the spin-dependent 
cross section in tDY. 

The factorization formula (2-3) provides not only a realization of the parton model based 
on QCD but also a systematic framework for calculating the QCD corrections beyond the 
parton model. As mentioned above, the meaning of the factorization formula (2-3) is the 
separation of the short- and long-distance contributions contained in the process (2-1) at 
the factorization scale fx F . In formal field theoretical language, (2-1) can be represented 
by Feynman diagrams, assuming that the hadrons h^ 2 can be described in terms of virtual 
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partonic states, as if they were a beam of quarks and gluons. Then, in a corresponding 
generic diagram, the long-distance contributions may give rise to infrared (IR) divergence, 
but, considering the set of diagrams contributing at the same order of accuracy, those IR- 
divergent contributions can be completely factorized into universal hadron matrix elements 
[transversity distributions, see (AT)], or cancel out, and thus the partonic cross section (2-5) 
possesses purely short-distance contributions containing all the dependence on the relevant 
hard scale Q. 1 ^ Therefore, QCD perturbation theory can provide a prediction for (2-5) in a 
systematic way, while the transversity distributions as the universal hadron matrix elements 
can only be determined by taking into account nonperturbative effects. 

A traditional approach for the separation of short- and long-distance contributions is the 
operator product expansion (OPE), and in fact, the OPE has been successfully applied to, 
e.g., deep inelastic lepton-hadron scattering. But, in the case of DY processes, as well as 
various hard processes in which two or more hadrons participate, the conventional OPE in 
terms of local operators is useless for this purpose, and the factorization theorems associated 
with formulae like (2-3) provide a generalization of the OPE suitable for those processes. 1 -* 
The partonic cross section (2-5) in the factorization formula (2-3) plays the role of the 
Wilson coefficient functions in the conventional OPE, and (2-5) can be determined by a 
"matching" procedure similar to that used to obtain the Wilson coefficient functions in the 
OPE: To obtain the prediction for (2-5), we apply the factorization formula (2-3) to the 
parton- level process (2-2), employing an appropriate IR regularization scheme. In this case, 
(2-4) can be calculated explicitly as the "quark distribution inside a quark," 5q hl=q (x 1 , /j? f ) = 
S(xi — 1) + • • ■ and 5qh 2 =q(x2, fJ>%) = o~(x2 — 1) + • • • , where the ellipses represent the terms 
of 0(a s ) or higher, depending on the specific IR regulator [see (A-3) in Appendix A]. The 
IR-regulator-dependent contributions on the RHS of (2-3) actually coincide with the IR- 
regulator-dependent contributions appearing on the LHS for the parton-level cross section, 
and thus a comparison of both sides of (2-3), order-by-order in a s , yields A F do-(p 1 ,p 2 ; Q 2 ; n F ) 
of (2-5) as an "IR safe" power series in a s . Through this matching procedure, we confirm that 
the partonic cross section, ATda(pi,p2~,Q 2 ; fi F ), indeed represents a purely short-distance 
component in the cross section (2-3). In our calculation presented below, we employ the 
dimensional regularization scheme to regularize IR and UV divergences. 

2.2. Momentum and spin for parton-level processes in D dimensions 

We now perform the above-mentioned matching at the one-loop level and derive A T da(pi,p 2 ] Q 2 ; [x 2 F ) 
in (2-5) to 0(ct s ). For this purpose, we calculate the LHS of (2-3) for the parton-level process 
(2-2), which we write as A T da qq . At the one-loop level, the diagrams to be calculated are 
given in Fig. 1. We calculate all these contributions in D = 4 — 2e dimensions to regulate 
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the IR and UV divergences. For these parton-level processes, we define the invariants as 



s = (pi + P2? , Q 2 = (h + k 2 f = zs , t = (pi - hf , u = (pa - h) 2 , (2-6) 

in terms of the momenta assigned in (2-2). In the case of one-gluon emission corresponding 
to the lower diagrams in Fig. 1, i.e., (2-2) with X = g(k), we also define 



t = ( Pl - k) 2 , u = - A;) 2 . 



(2-7) 



The momentum of a massless particle, p = {pi,P2, k±, k 2 , k}, in D dimensions is generically 
expressed as 



p" = (P° , P) = (|p| , P), P = \P\H D ~ 2) sm9 p , cos0 p ), 



(2-8) 



where rip D ^ denotes a (D — 2)-dimensional unit vector which is parametrized by D — 3 
angular coordinates, 4> P D ~ 3 \ (pp '^, . . ., as 



n 



( D - 2 ) = (cos0l D - 3 ), nl^sin 



-3b 



= (cos0f- 3 ), sin^f^cos^f- 4 ), nf- 4 )sin0( D - 3 )sin0f- 4 )) 
. . . , 

such that in 4 dimensions, defining (f)^ = p , we have 

n£ 2) = (cos0 p , sin0 p ). 



(2-9) 



(2-10) 





P2 




wftrtrtrcffffcrffip 
k 



wfffowGWtrtP 




Fig. 1. Parton-level processes contributing to the spin-dependent cross section in tDY to 0(a s ). 
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When we calculate the LHS of (2-3) for the parton- level process (2-2), we always choose a 
frame in which the initial quark and antiquark in (2-2) possess momenta p x and p 2 that are 
collinear along the z-axis, with 9 Pl = and 8 P2 = ir, and possess transverse spins, satisfying 
Pi • si — 0, p 2 ■ s 2 = 0, and s\ 2 = — 1, that are expressed by 

< 2 = (0, s[ D ~ 2 \ 0), (2-11) 

with the (D — 2)-dimensional unit vectors s[ D 2 ~ 2 \ which reduce in 4 dimensions to 

= (cos 0i, sin0i), = (cos 02, sin0 2 ). (2-12) 

2.3. Parton-level amplitude squared 

The contribution of each diagram in Fig. 1 to the spin-dependent cross section Arda^, 
the LHS of (2-3), is generically given by 

^^t|M| 2 ^, (2-13) 

where d<P denotes the differential element of the n-body phase space for the corresponding 
2 — > n process (n = 2,3), 1/N C is the color-averaging factor for the incoming quark (an- 
tiquark), and the spin-dependent amplitude squared, Ar|M| 2 , is defined in terms of the 
relevant 2 — > n Feynman amplitude M(si,s 2 ), associated with the transverse spins s\ and 
s 2 for the colliding quark and antiquark in (2-2), as 

A T \M\ 2 = l - [|M( Sl ,s 2 )| 2 - |M( Sl ,-s 2 )| 2 ] , (2-14) 

where the sum over the spins of the final-state leptons is implicit. We calculate (2-13) and 
(2-14) in the dimensional regularization. We rewrite the QCD coupling constant as g — > gfi e 
with the mass scale \i. This enables us to keep the redefined coupling constant dimensionless 
in D = 4 — 2e dimensions.*) 

In the spin-dependent amplitude squared of (2-14), we encounter traces of the gamma 
matrices involving 75. We employ the naive anticommuting-75 scheme, {75, 7 M } = (// = 
0,1,... D — 1), which is the usual prescription in the transverse-spin channel: For the 
transverse-spin case, it appears that these traces involve only even numbers of 75. Therefore 
the matrices 75 eventually disappear due to 7! = 1, and we do not anticipate any inconsis- 
tencies related to 75. This suggests that the naive anticommuting-75 prescription will work 
(see Ref. 19)). 

*) In principle, one can rewrite the QED coupling constant similarly. This results in the multiplication of 
the Born cross section and the higher-order cross section by a common factor of /i 4e , and it allows us to keep 
the dimensions of the cross section unchanged. But this is irrelevant for the present purpose of matching to 
obtain the partonic cross section (2-5). 
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The upper two diagrams in Fig. 1 represent the 2 — > 2 processes up to 0(a s ). The 
corresponding tree-plus- virtual contributions to (2-14) can be calculated straightforwardly, 
and they read 



A T \M\ 2 T+Y = 4N C 



( i.Tiae q 



2 s (ki • si)(k 1 ■ s 2 ) +tu(si- s 2 ) + e— (si ■ s 2 ) 



x 



a. 



1 + —C F 



7T 



r l-e 



(2-15) 



where CV = (iV 2 — l)/(2N c ) [also see (2-6)]. Equation (2-15) includes the contribution 
due to the wave-function renormalization factors for the incoming quark and antiquark legs, 
which cancels the quark-photon vertex renormalization constant, as guaranteed by the Ward 
identity Therefore, (2-15) is UV finite, and it is also gauge invariant in a gauge-invariant 
regularization scheme. Note that the ratio of the 0{a s ) term to the 0(1) term in (2-15) is 
the same as the corresponding ratio in the unpolarized case. 40 ) 

The lower two diagrams in Fig. 1 represent the 2 — > 3 processes up to 0(a s ). The relevant 
polarization sum for the final-state gluon, V a/3 (k) = J2x e (A)(^) e (A)(^)' wnere e (\)(^) * s ^ ne 
polarization vector for the real gluon with momentum k and the physical polarization A, 
can be replaced as V a/3 (k) — > — g af3 , using the Ward identity, in the present case of one real 
external massless gluon. 43 ^ We write the corresponding result for the one-gluon-radiation 
contributions as 



A T \M\l = N C C F 



Q 2 



(2-16) 



where AtTZ is formally given by a second-order polynomial in e, after working out the traces 
of the Dirac matrices. Actually, the terms proportional to e 2 in A^R, do not contribute to 
the cross section in the limit of e — > 0, and therefore we drop them. Similarly, we drop some 
of the terms proportional to e, which clearly do not contribute to the cross section, and we 
thereby obtain 

16e 

A T K = — [s (zt + u + Q 2 )(h ■ s 2 )(k ■s 1 )-tu(k- s r ){k ■ s 2 )] 

6 

86 

+-r [2 (s - Q 2 ) (h ■ si)(fci • s 2 ) + u (t + u + Q 2 ) (si • s 2 )] 

t 

16e 

+— \s(zu + t + Q 2 )(h ■ Sl )(k ■ s 2 )-tu(k- Si)(k ■ s 2 )l 
u z 

8e 

+— [2 (s - Q 2 ) (h ■ si)(fci • s 2 ) + t (u + t + Q 2 ) (si • s 2 )] 



16Q 2 



g2 s 

2 s (ki ■ s 1 )(k 1 ■ s 2 ) + iu(si • s 2 ) + e-^-(si • s 2 ) 
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80 2 

[2 s {(h ■ s 2 )(k ■ si) + (fci • si)(A; • s 2 ) + (A; • si)(fc • s 2 )} 

+ {Q 2 (s + t + M) -tt-MM + tM}(si ■ S 2 )] 

8e 

+— s (i + u + Q 2 ) (s + t + u) (si • s 2 ) 
iw 

+i [i (s + * + u) + (s + 1) (t + u + Q 2 )] (si • s 2 ) 
4e 

+— [u(s + t + «) + (s + u) (w + t + Q 2 )l (si-s 2 ), (2-17) 

■u 

using the variables in (2-6) and (2-7). Note that the terms in (2-17) may receive additional 
factors of 1/e as IR-divergent poles through the integration over the phase space in (2-13). 

§3. NLO matching for the Q T -integrated (total) cross section 

First, we consider the case in which the transverse momentum, Qt, of the final dilepton 
is unobserved (integrated) in the parton-level process (2-2). We calculate the corresponding 
cross section (2-3) using (2-13)-(2-17) and determine the (^T-integrated partonic cross section, 
A T da( y p 1 ,p 2 ;Q 2 ;fi 2 F ), to 0(a s ). 

3.1. Tree-plus-virtual contributions 

The phase space for an outgoing massless particle with momentum p is given by [see (2-8) 
and (2-9)] 

(2tt)^- 1 2/ 2(2tt) d - 1 M M p ' 

where the differential element for the (D — 2)-dimensional angular integration is expressed 
as 

dOf-V = (1 - cos 2 # p )- e dcos# p dflf- 2) . (3-1) 
The two-particle phase space for the 2 — > 2 processes of Fig. 1 is given by 

* - &Pk(S^k (2I)D * ,D) (pi +K - ki ~ k2h (3 ' 2) 

and we obtain for the tree-plus- virtual contributions to (2-13) the expression 

using the invariants of (2-6). Inserting the expression for Z\t|M|t +v , (2-15), and using the 
explicit forms (2-8)-(2-ll) for the lepton's momentum k\ and the quark's transverse spin Si j2 , 
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we obtain 



A T da™ 



T+V 



/ 4-7rae r 



(D- 



dQ 2 dQ£ 



!) 16tt 2 AT c V Q 2 

1 -2/. 6 

-a* sin 9 kl + - (si ■ s 2 ) 
4 2 

a 



16tt 2 J 



S(l-z) 



X 



X 



r(i-c) 



1 3 7T ... . 

-?-2?- 4 + T + c,(£) 



where 



with 



a<2> = 2 (n 



«1 • «2 = - S 



(D-2) (D-2) 



(3-3) 



(3-4) 



(3-5) 



From (3-1), the corresponding "azimuthal" angular distribution of the lepton is obtained by 
integrating (3-3) over the "scattering" angle 9^, with the measure (1 — cos 2 flfcj -6 dcosO^, 
as 



A T da qq 



T+V 



( 



47rae r 



6(1- z) 



dQHQ { ^ 2) i^ 2 n c V Q 2 J VievrV 

x [2 1 - 2e 5(2 -e, 2-e) c^ + e2- 2e 5(1-6, 1 - e) ( Sl • s 2 )] 

■ 2 ^ e 1 / 1 3 7T 2 ,^ , 



X 



(3-6) 



where B(x, y) = r(x)r(y)/ T(x + y). 

3.2. i?ea/ g/non emission 

For the 2^3 processes depicted in Fig. 1, the real-gluon emission contributions, the 
corresponding amplitude squared seems rather complicated [see (2-17)]. Fortunately, mirac- 
ulous cancellations occur among the first four lines and among the last three lines of (2-17): 21 ) 
Because the corresponding terms are proportional to e, these terms survive only in the con- 
figuration obtained in the t,u — > limit, in which the divergences associated with the 
collinear (k^ 1 oc pi, k^ oc p%) or soft {k^ 0) gluon are produced [see (2-7)]. An explicit 
calculation shows that the collinear configuration is indeed relevant in the present case, and 
in this configuration an additional 1/e factor multiplying these terms appears. However, 
for such collinear configurations, it turns out that the first (third) line cancels the second 
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(fourth) line in (2-17), and also, the last three lines cancel. Therefore, (2-17) can be reduced 
to the very simple form 

A T TZ = A T 1ls + A T n NS , (3-7) 

with 



A T n s 

A T 1Z NS 



16Q 2 
tu 
SQ 2 



Q 2 s 

2s (h ■ si)(*4 • s 2 ) + tu(s 1 ■ s 2 ) +e—^-(si ■ s 2 ) 



2s {(h ■ s 2 )(k ■ si) + (fci ■ si)(/c • s 2 ) + (A; • s 1 )(k ■ s 2 )} 



tu 

+{Q 2 (s + t + u) -ti-uu + iu} (si • s 2 ) 



(3-8) 



We have divided the RHS of (3-7) into two terms; as will be discussed below, only AtTZs, 
given in (3-8), receives the 1/e singularity through the integration over the phase space in 
(2-13), while AxTZns yields the result finite in the limit e — > 0.*- 1 

The three-particle phase space for the 2 — > 3 processes depicted in Fig. 1 is given by 

* = (ffk (2^*8 (2^° (2 *» D ^ *")■ <^> 

To use this in (2-13), it is convenient to employ the qq CM frame and first determine the phase 
space factor corresponding to the differential elements, dQ 2 dQ^ ^ df2^ x \ integrating 
over the other degrees of freedom [see (2-6)]. This results in the replacement 
1 Q2(i- 2e ) Ik] 1 - 2 * 



d&i 



where 



d n 2 d<l n 1 d<> " 

8^(2vr)^ 2 ™ (QO + fc.fcVlfcil) 2 -*^ dJ4 



Q° = ^(l + *), \ k \ = £(l-z), 



(3-10) 



and the magnitude of the lepton momentum, = (|fei| , fei), is fixed by the relevant kine- 
matical constraints as 

Q 2 l __ Q 2 l 



2 QO + fc.fcx/lfcil ®i 



(3-11) 



with [see (2-8)] 



0/ = 1 + z + (1 - s) (cos fl fc cos fcl + sin fl fc sin fcl njj? ~ 2) • njf~ 2) ) , (3-12) 
which is positive definite. The invariants in (3-8) are now given by [see (2-6), (2-7)] 



* = --(l-*)(l-COS0 fc ) 

t = — \fs \ki\ (1 — cos^J, 



u = - - (1 - z) (1 + cos6» fe ), 

U = — | fci| (1 + COS flfcj. 



(3-13) 



*) Note that we have k ■ Si,2 ~ ► and Q 2 (s + 1 + u) — ti — uu + iu — > in the collinear limit, that is, 
fc^ || or fc^ || p%. 
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Substituting the above results, the cross section (2-13) for one-gluon emission in the qq CM 
frame reads 

A T dag 



(D-l) 



where we have 



G(e) 



dQ 2 df2 y ki 
Q 2 

' 6A7r 2 N r s 2 



G{e) 



\l-2e 



A V ( X ~ Z "> jnl - 1 ) 



(3-14) 



/ 4nae q x 2 



Q 2 



Q 2 



16tt s 



a 

TV 



s c F 



~Q 2 



2\ e 



7T C 
47r' 



With (3-11), (3-13), (3-4), and (3-5), the singular part of A T K, i.e., A T TZ S given in (3-8), 
reads 

16 Q 2 



A T K S 



tu 
64 Q 2 



-—«<?> sin kl + e — - 



si ■ s 2 



1 



1 -z) 2 sin 2 ^ L^/ 



g 2 s 



a<2> sin fc + e — — S\ ■ s 2 



Using (3-1) for the angular integration of the gluon in (3-14), and changing the integration 
variable cos^ to z\ = (1 + cos^) / 2, we obtain the contribution to (3-14) from the singular 
part of one-gluon emission as 



A T daf {s) 



dQ 2 dQ K k \ 



^ = G(e) 



32 Q 2 



-2t 



(D-l) 



s 2 (i- z y+ 2 * 



J dz x J dQf- 2) 



x 



Q 



— a$ sin kl + e ■ 



,(l-^i) 1+< 

Q 2 s 



+ 



(l-*i)- 



T 2-2e 



Si ■ S 2 



(3-15) 



Because @)i depends on the angles of both the gluon and lepton as expressed in (3-12), it is 
difficult to work out the integration of (3-15) in D = 4 — 2e dimensions. However, fortunately 
we can carry out the integration by noting the well-known formula 



1, . 1 



ln(l - z) 
l-z 



(3-16) 



and the similar formula obtained through the replacement z — > 1 — z in (3-16); these formulae 
allow us to make manifest the singularities appearing in the limit e — > 0. Here the "+" 
distributions are defined, as usual, as 



dzf(z) 



(l-z) 
\n(l-z 



+ 



1 — z 



(3-17) 



for a smooth test function f(z). Note that the 1/e poles proportional to 5(1 — zi) and 
when employing the formula (3-16), are associated with the radiation of the collinear gluon 
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from an initial particle as z\ — (1 + cos#fc) / 2 = 1 or 0, while the 1/e pole proportional to 
8(1 — z) is associated with the radiation of the soft gluon at the threshold for the reaction, 
z = Q 2 /s = 1. After inserting (3-16) into (3-15), we need to carry out the integration in 
arbitrary D = 4 — 2e dimensions only for the terms proportional to 1/e 2 and 1/e, associated 
with the following limiting cases for 

® I (z^l) = 2, (3-18) 
9i(z x -> 1) = l + z + (l-z) cos6 kl , (3-19) 
^j( Zl _>o) = l + z-(l-z) cos8 kl . (3-20) 

Then, the corresponding angular integrals in (3-15) can be performed exactly. The integra- 
tion of the remaining terms in (3-15), which are finite as e — > 0, can be performed in D = 4 
dimensions, using [see (2-10)] 

9i(D -> 4) = 1 +z+ (1 - z) (cos# fc cos9 kl + sin^ fc sin# fcl cos (0 fc - <f> kl )) . (3-21) 

The corresponding angular integrals can be done straightforwardly, by using the formulae in 
Appendix A of Ref. 16). Therefore, the phase space integrals in (3-15) are tractable exactly 
up to corrections that vanish as e — > 0. 

We write the contribution to (3-14) from the nonsingular part of A T 1Z, i.e., A T 1Z NS 
of (3-8), as A T da^ NS ^/ dQ 2 df2^~ 1 \ Similarly to the finite terms in (3-15) discussed just 
above, the relevant angular integrals in A T da^ N ^ / 'dQ 2 'dQ^~^ can be calculated in D = 4 
dimensions, where AtTZns reduces to 



A T 1Z NS 

32 Q 2 



[1 - z) 2 sin 2 Qi 



'1 - z) 2 

- — r sin 2 6 k cos(0 fc - 0i) cos(0 fc - 2 ) 



+Z(1 - Z) Sm 6 ' fcl {cos(0 fe - 0i) COs(0 fcl - 02) + (1 <"> 2)} 

. (l-^) 2 .2/, n . x 1 / n .cos^fccos^l 
2 + ^ sm # fe - z(l + z) — - z(l - z) — I cos(0i 



(3-22) 



with 12 parameterizing the spin vectors Si i2 in 4 dimensions as in (2-12). Using (3-22), 
(3-21) and the formulae in Appendix A of Ref. 16), the calculation is straightforward, and 
gives At da^, NS J l dQ 2 dQ ( k ® l \ up to corrections that vanish as e — > 0. 

Combining the above results, we obtain At da^/ dQ 2 df2 ( k ® ^ = 
A T daf {s) /dQ 2 dQ ( k D i ~ 1) + A T daf {NS) /dQ 2 dQ [ k D i ' 1) . Then, using (3-1) with p -> ^ and inte- 
grating over # fcl , we finally derive the "azimuthal" angular distribution of one of the leptons 
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for one-gluon radiation: 

16 vr 1 - 6 



A T da™ 



dQHrt»- 2) r(i - e) 



G{e) Q 2 



-J{l-z)-- 



e (l-z) + 



x [2 1_2e B(2-e, 2-e) a* + e2- 2e B(l - e , 1 - e)( Sl ■ s 2 )] 
+ 16 7t G(0) Q 2 | cos(20 fcl - 0x - 2 ) 



x 



7T 2 r/i , /lnl-z \ l-z 3 In 2 z lnz 

6(1- z) + 2 — : H 

12 v ; V l-z J , z 2 l-z l-z 



(3-23) 

Apparently, this final form is valid in any frame in which the momenta p 1 and p 2 of the 
initial quark and antiquark in (2-2) are collinear along the z-axis. 

3.3. Partonic cross section to 0(a s ) in the MS factorization scheme 

With the tree-plus- virtual contribution (3-6) added to (3-23), the double poles in e cancel 
out, while the single-pole terms remain. The result gives the azimuthal angular distribution 

to 0(a s ) as 



A T da qq ~ 



A T da™ 



T+V 



A T da™ 

,(D-2) 



+ 



dQ 2 dQ i f i ~ 2) dQ 2 dd®- 2) dQ 2 dfl£~ z> dQ 2 df2£- z> dQ 2 dQ k \ 



A T da™ 

,(D-2) 



+ 



A T da™ +R 



(D-2) ' 



(3-24) 



for the parton-level process (2-2) in the dimensional regularization. Here, in the second 
equality, we have rearranged A T da qq / dQ 2 df2 k ® ^ into a sum of the O{a s ) and 0{a\) con- 
tributions. The remaining single-pole terms in At da^ +K / dQ 2 dQ ( k ®~ 2 ' S> represent the mass 
singularity associated with the emission of the collinear gluon. Indeed, these single-pole 
terms are completely absorbed into the parton distributions as the corresponding mass sin- 
gularities: We substitute the result (3-24) into the LHS of the factorization formula (2-3) 
applied to the parton-level process (2-2) and (A-3) in Appendix A into 8H{x\,X2] fi 2 F ) in the 
RHS as the transversity distributions in the MS scheme for the incoming quark and anti- 
quark [see (2-4)]. Matching both sides of the resulting formula, we determine the partonic 
cross section A T da/dQ 2 df2 k ^~ 2 ^ on the RHS order-by-order in a s . The matching at LO 
shows that the O(a s ) term of A T da j dQ 2 dQ^ 2 ^ reads 



A T da^( Pl ,p 2 ;Q 2 ;fi 2 F ) 



A T da™ 



dQ 2 dQ ki 



(D-2) 



n dQ 2 dn k \ 



(D-2) 



16tt 2 A" c V Q 2 

x [2 1 ~ 2e B(2 -e, 2-e) c^ + e2- 2e B(l-e, 1 - e) (s 1 ■ s 2 )] , (3-25) 
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exactly for arbitrary e. Apparently the RHS of (3-25) is free from the 1/e singularities and 
does not depend on fip- Similarly, the matching at the NLO level yields the 0(a l s ) term of 
A T da/dQ 2 dn { °- 2) as 



A T d&^( Pl ,p 2 ;Q 2 ;^ F ) _ 1 A T da™ +R 



+ Yf^ j dxA T P qq {x) 



dQ 2 dQ { °- 2) dQ 2 d(l£ 

I At da^(xp 1 ,p 2 ;Q 2 ;fi 2 F ) A T da (0 \p 1 ,xp 2 ;Q 2 ; /j, 2 f) i 
1 dQ 2 drt°- 2) dQ 2 drtt 2) ' 



for the MS factorization of the mass singularities at the scale fj,F, where 1/e = 1/e — ^ E + 
ln(47r), with 7^ the Euler constant, and A T P qq {x) is the LO Dokshitzer-Gribov-Lipatov- 
Altarelli-Parisi (DGLAP) splitting function for the transversity distributions, given by (A-4). 
The terms proportional to 1/e in (3-26) are generated by the 0(a s ) term of the parton 
distributions (A-l), combined with the O(a s ) terms of the partonic cross section, (3-25), 
and cancel the 1/e poles in the first term on the RHS of (3-26). As a result, (3-26) is also 
finite as e — > 0. 

Now, the mass singularities for the parton-level process (2-2) have been completely fac- 
torized into the relevant parton distributions, and the e — > limit of the sum of (3-25) and 
(3-26) gives the final result for the partonic cross section (2-5) to 0(a s ) in 4 dimensions, 

A T da(p l ,p 2 ;Q 2 ;^ 2 F ) a 2 

x(i<i-, ) + £ [2 ,„J + Cf { & (!^) 



„zhi 2 z z\nz w „ . f2n 2 
- 6 T -j-A^— + A(l-z)+( — -8)5(1 -z) 



+ 

, (3-27) 



where 0^ and <pi y2 are defined as in (2-10) and (2-12), [i F is the MS factorization scale, and 
we have used the fact [see (3-4)] 

a<p(D ->• 4) = cos (20 fel - 0i - 4> 2 ). 

The result (3-27), along with (A-4), coincides with that obtained in previous works employing 
the massive gluon scheme 16 -* and the dimensional reduction scheme, 17 - 1 via the scheme trans- 
formation relation 19 ) (see also Ref. 20)). Substituting (3-27), and also the NLO transversity 
distributions for the hadrons h± t2 in the MS scheme, into the RHS of (2-3), we obtain the 
NLO QCD prediction for the spin-dependent cross section of tDY (2-1) in hadron-hadron 
collisions, A T da/dQ 2 d(p kl , as the mass (Q) and the azimuthal angular (0&J distribution 
associated with the observed lepton pair. 
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§4. The Q T differential cross section to 0(a s ) 

In this section, we extend the calculation of the last section to the case in which the 
final state in (2-1) is observed in more detail. We calculate the spin- dependent cross section 
for tDY, which is differential also in the transverse momentum Qt and rapidity y of the 
produced lepton pair, 

^ (4-1) 

dQ 2 dQ 2 T dyd^ 1 ; 

taking into account the QCD mechanism up to 0(a s ) in perturbation theory. The rapidity 
of the lepton pair is defined by 

in the hadron-hadron CM system for (2-1), in which we work in this section. Here, Q° and 
Q 3 are the components of the dilepton's momentum, = (Q°,Q T ,Q 3 ), with Q 3 denoting 
the component along the direction of the colliding beam, while Qt is defined as Qt = \Qt\- 
Now, in the factorization formula (2-3) for (4-1), the corresponding partonic cross section 
on the RHS depends also on Q T and y, as A T da(p 1 ,p 2 ; Q\, Q 2 , y, n%) /dQ 2 dQ T dy d<f> kl ■ 
Similarly to the last section, this partonic cross section can be obtained by matching the 
LHS and RHS of the factorization formula (2-3) in the case of the parton-level process 
(2-2) with Q T and y observed. The LHS including the 0(a s ) corrections is represented 
by the diagrams in Fig. 1, which we calculate using dimensional regularization and obtain 
A T da qq / dQ 2 dQ\ dy dQ^ 2 ^ = A T da^ +v+R / dQ 2 dQ T dy dft^ 2 ^ as the tree-plus- virtual con- 
tribution plus the real emission [compare with (3-24)].*) As demonstrated in the last section, 
the corresponding order-by-order matching accomplishes the factorization of mass singu- 
larities in terms of 5H(xi, x 2 ) of (2-4). Because 5H(xi, x 2 ) is universal (process indepen- 
dent), the factorization of the mass singularities can be expressed by the relation between 
A T da q ydQ 2 dQ 2 T dydn { °~ 2) and A T da / 'dQ 2 dQ\ dy dff k D ~ 2) ', which is analogous to (3-25) 
and (3-26): 



A T da(p u p 2 ;Q T ,Q 2 ,y;fx 2 F ) 1 



dQ 2 dQ 2 T dydQ ( k D i ' 2) ~~ e 2 q 



A T da qq ~ 



T+V+R 



dQ 2 dQ T dydn K k \ 



(D-2) 



A T daf 



+ 'CT 



dQ 2 dQ\dydnf- 2) 



_ A T rffrr+v+R+CT(pi,P2; Q T ,Q 2 , y; /4) ^ 4 . 3 j 



dQ 2 dQ 2 T dydni\ 



(D-2) 



*) Here and below, Axdcy /dQ 2 dQ\dydfl k denotes the contribution of the virtual-correction di- 
agram in Fig. 1, combined with the contribution due to the wave-function renormalization factors for the 
incoming quark and antiquark legs [see the discussion below (2-15)]. 
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where the second term on the RHS of the first line plays the role of the "counter-term" to 
cancel the mass singularities in the "V+R" contributions contained in the first term, and it 
reads, in the MS factorization scheme, 



At da, 



qq 

CT 



dQ 2 dQ 2 T dydQ^ 



(D-2) 2e 



\_Os_ r 

le 7T J 



dxA T P qq {x) 



x 



A T da™ 



dQ 2 dQ T dydn£ 



(D-2) 



+ 



A T da™ 



T 



dQ 2 dQ T dydn£ 



(D-2) 



(4-4) 



with the LO splitting function A T P qq (x) of (A-4). From (2-13), (2-15), and (3-2), it is 
straightforward to show 



A T da™ 



T 



S(Qt)S y 



2 



A T da™ 



(4-5) 



dQ 2 dQ 2 T dydd^~ 2) -< y 2 PiJ dQ 2 dfi)^ 

exactly in arbitrary D = 4 — 2e dimensions, where At da^ / dQ 2 dQ^ ^ on the RHS is 
given by (3-25). At the tree level, apparently, it is impossible for the outgoing dilepton to 
have nonzero transverse momentum, and the rapidity is also determined completely by the 
kinematics of the initial state given by (2-2). 

It is not difficult to identify the 1/e poles in the "V+R" contributions through a careful 
treatment of the relevant phase space integration in D dimensions, as done in (3-16)-(3-21) 
of the last section, and in so doing to verify that those pole contributions are indeed canceled 
out by the "CT" contribution in the partonic cross section (4-3). The result would yield the 
partonic cross section, which is finite as D — > 4 and is expressed in terms of the relevant 
partonic variables. However, rather than performing such manipulations directly at the 
partonic level, it is actually more convenient to perform the corresponding manipulations in 
the factorization formula (2-3) for the hadron-hadron collisions, with (4-3)-(4-5) substituted 
into the RHS as it is and their partonic variables re-expressed in terms of the relevant 
hadronic variables of (2-1). Because the finite value of the dilepton's transverse momentum 
Qt is a consequence of the recoil from the gluon radiation, the IR behavior associated with 
the real gluon emission is controlled by Qt- Thus, rearranging the convolution integrals 
in the factorization formula (2-3) to make explicit its behavior as Qt — > 0, we confirm 
the cancellation of the 1/e poles and simultaneously get the cross section formula for the 
hadron-hadron collisions, which is organized according to the Qt dependence; such a form 
of the cross section is particularly suitable for the calculation of the dilepton Qt spectrum 
in hadron-hadron collisions, and it is also even desirable when we attempt to include the 
soft gluon resummation effects relevant in the small-Qr region in the next section. 
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We denote the total CM energy in (2-1) as y/S, and thus we have Pi = |v / S'(l,0, 1), 
P2 — |v / 5'(l, 0,-1), where denotes the D — 2 dimensional null vector. We decompose the 
cross section (2-3) into four pieces, corresponding to the contributions on the RHS of (4-3), 
as 

A T da = Z\ T d<7 T +v+R+CT, (4-6) 

and introduce the following useful hadronic variables according to the treatment of the 
unpolarized DY process in Ref. 31): 



S 



T — — — , %-\ — \j T 6^ , Xn 



Xi = (-5-J e • 12 = ~vr = l^H e ■ (4 ' 7) 

Here, the definition (4-2) of y is used, and the DY scaling variables, as usual. First, 

using (2-13), (2-16) and (3-9), the partonic cross section for the one-gluon radiation of (4-3) 
can be written for pi = XiPi as 

A T d& R ( y x 1 P 1 ,x 2 P 2 ) = A T da^\ pi=XiPi 

(4-8) 

where, for simplicity, we have denoted the factorization scale as fi and have suppressed 
some arguments of the partonic cross section Axda^- We have also introduced the following 
shorthand notation: 



7P 



^ = TITT = Q° - cos ^i " Qt *n9 kl n%- 2) ■ n^- 2 \ (4-9) 
l^i I 

with Q ' 3 = (Q 2 + Qr) 1/2 (e y ± e- y )/2 from (4-2). The quantity @ T is generated in (4-8) as 
a phase space factor from while the delta function 5 {{x\P\ + X2P2 — Q) 2 ) comes from 
the on-shell condition for the final-state gluon and is given by S(s + t + u — Q 2 ) in terms 
of the partonic variables in (2-6) and (2-7). Using the hadronic variables (4-7), this delta 
function becomes 

5 ((xiPi + x 2 P 2 - Qf) = 1 5{x lX2 - Xl x+ - x 2 x+ + r). (4-10) 
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It is important to notice that, in general, a convolution integral like (2-3) with the delta 
function (4-10) can be split into the two symmetric integrals 31 ^ , *) 



dx\dx 2 f{x\ , x 2 ) S(x\X 2 — X\X 2 — x 2 x\ + r) 



-/ 

J .Ft 



fjXj , X * 

axi 

y/T^eV X\ — X 



1 Kn 



dx' 



^e-y 



+ 1 

x 2 - xj 



where**) 



3j O 



tJU ]^ t,L 2 t,t 1 tX. 2 



Qt 



+ \/T 



In A T TZ of (4-8), we first consider the singular part of (3-8), which now reads 

A T lls Q 2 



16 Q 2 



Qt 



Q 2 l 

a$ sin 2 6 kl + e- (si • s 2 ) 



and this yields the following contribution to (4-8): 



A T da R t S )(xiPi , X2P2) 



'R(g) 

dQ 2 dQ T dy 



= G 



Hi 



1 



' ' (Q 2 n) 1 + e X X X 2 

^b[x\x 2 — x\x\ — x 2 x\ + t) df2g 



-M T (e,Q 2 T ) 



^df2^\ 



(4-11) 



(4-12) 



(4-13) 



(4-14) 



where 



l Q 2 



M T (e, Q| ) = — a<P sin 2 + . ^ 



1 ^ ^ 
e-(si • s 2 ) . 



(4-15) 

Inserting (4-14) into (2-3) and using (4-11), we obtain the contribution to A T da R in (4-6), 
associated with the singular part (4-13), 



A T da R{S ) _ Gli 



1 



dQ 2 dQ T dy 



dXl XTTf *\ T 

■ on [x\ , x 2 ) 



e y X\ X l 



+ 



f 1 dx 

Jy/T^e-y X 2 — 



2 — 5H(x\ , x 2 ) -J— 



X\X 2 

M T (e,Q 2 T ) dng-^dng-v, 



(4-16) 



where we have suppressed the /i 2 dependence of 8H(xi , x 2 ). We follow the procedure used 
for the unpolarized DY process in Ref. 31), in order to isolate the singular terms in the limit 



*) When ^/r+e y > 1 (yfr+e~ v > 1), the first (second) term in the RHS of (4-11) vanishes, because of the 
support property of the parton distributions contained in f(x\,X2) [see (2-3)]. 
**) Note that we have xf 2 = ^ff^e ±y = x\ 2 = %\ 2> when Q\, = 0. 
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Q\ — > of the above equation, corresponding to the radiation of the collinear gluon: 

A T da R ( S ) 



dQ 2 dQ\,dy 
= G H (e 



+ 



(Qt) 1 ' 



e y X\ — X 1 



8H(xi , x\ 



SH(4 , x° 2 ) 



x\x 2 



dxi 



/t + e~» x 2 x 2 

+ 5H(x° 1 , xl) In 



8H{x\ , x 2 ) 



T 



X-t X2 



- 5H{x\ , xl) 



Jh J\ .r,)ll r 2 )S 
Qt 



(D-2) , n (D-l) 



M T (e, Q£)df% ; <W2£ 



(4-17) 



We can further isolate the poles in Q\ by using the identities 



Q 



2c 



(<5t) 



2 U+e 



2c 



ln^ 

(Qt) 1+€ Qt 



77^r-^(Q 2 T ) + 0(e), 
Wt)+ e 

HQ 2 /Q 2 i 



±6(Q T ) + 0(e), 



(4-18) 



Qt J + e ^ 

where the "+" distributions that regulate the singularity at Q\ = are defined such that 32 ) 

'\n n (Q 2 /Q 2 T )\ I 



n + 



-\n n +\Q 2 / V 2 T ). 



(4-19) 



Then, analogously to the situation discussed with regard to (3-18) - (3-21) for the case of 
the Qr-mtegrated cross section, we find that only the two limiting cases for M T (e, Q 2 -), 

M T (e = 0,Q 2 T ^0) and M T (e ^ , Q 2 T = 0), 

are required in order to evaluate the (D — 3)-dimensional angular integration dQ^ 2 -* in 
(4-17), up to corrections that vanish as e — > 0. In these two cases, the integration over 
the scattering angle 8 kl of the lepton can be also performed easily. After straightforward 
calculations, we obtain the "azimuthal" angular distribution of the lepton for one-gluon 
radiation as [see (2-10), (2-12)] 



A T daji(s) 



2ir 



^- G H (0) cos {2<t> kl - 0i - fa) ( 1 + 



dQ 2 dQ 2 T dydf2^' 2) Q 



2Q\ 
Q 2 



x <J 5H{x\ , x° 2 " 



(Q 



1 ln (l-xt)(l-x+)S | f HQ 2 /Q 2 T ) \ 
t)+ Q 2 V Qt J 



dxi 



+ 



y^fev x l x l 

' dx 



+ 



5H{x\ , x 2 ) 



X\X 2 



- 5H{x\ , xl 



(Qt)+ J ^/T^e-v x 2 x t 



5H(x\ , x 2 ) 



x\x 2 



- 5H{x\ , x%) 
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2-7T 7T~ e 2~ 2e 

+ 7V2 Tn V [ 2 5(2 - e , 2 - e) ^ + e 5(1 - e , 1 - e)( Sl • s 2 )} 5(Q 2 T ) 



x <J 5H (x? , x^ 



1 _1 ln (l_-x?)(l-^ 



e 2 e 



o Xi — X" 



Q 2 



SH( Xl ,x%) — - 5H(x\ , x° 2 ) 



+ 



f 1 dx 2 

Jx9, x 2 — x 2 



5H(xl ,x 2 ) — - 5H(xl , x° 2 ) 
x 2 



(4-20) 



It is straightforward to show that the partonic cross section for the virtual correc- 
tion, A T da Y (p u p 2 )/dQ 2 dQ 2 T dydf2 { °~ 2) of (4-3), is given by (4-5) multiplied by (a s C F /e 2 w) 
(Anfi 2 /Q 2 Y[l/r(l - e)][-l/e 2 - 3/(2e) - 4 + tt 2 /2 + (9(e)], which is the factor appearing in 
the last parentheses of (3-6). Substituting this into (2-3), the virtual-correction contribution 
in (4-6) reads 



Arpda 



v 



dQ 2 dQ 2 T dydn { °- 2) 
2ir ir~ e 2~ 2t 

= 7V2 ~Fn V [ 2 S(2 - e , 2 - c) a, + c S(l - c , 1 - c)( Sl • s 2 )] 5(Q 2 T ) 



x <5#(x^ , x^) 



1 3 . 7T 2 v 

- ? - S -4 + y + 0( e) 



(4-21) 



Combining this result with (4-20) and using the identity 



JxO L — z J x o (1 - ZJ + 



we get 



A F da 



R(S)+V 



dQ 2 dQ 2 T dydf2®- 2 Q 



Q 2 



x ^ 5#(x" , x" 



(Q 



i ln (i-x+)(i-x+)g | / in(Q 2 /<&n 



(ixi 



(Qt)+ J yJT^eV x l — x i 



5H(xi , xJj) 



37 ]_ $2, 



5H(x° , x°) 



+ 



1 



(Ql 



dxo 



+ ^ G H (e 



T)+ J^/r^e-y X 2 x 2 

7i- e 2~ 2t 



5H(xl,x 2 )- <5#(x°,x°) 

Xl x 2 



[2B(2-e,2-e) a$ + e5(l-e, 1 - e)(si • s 2 )] <5(<2r) 
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X 



7T 



2G F e 



(4-22) 



with A^P^rr) of (A-4). We now derive the "CT" contribution of (4-6) in the MS scheme, 
using (4-3)-(4-5) in (2-3), and add the result to (4-22). We observe that the mass singularity 
poles in e arising in (4-22) are completely canceled by the 1/e poles of the CT contributions, 
and taking the e — > limit, we obtain in 4 dimensions 



A T da R ( S )+v+CT 2tt „ 

= — G (0) cos (20 fcl - 0i 



dQ 2 dQ^dyd4> kl Q 2 
1 



^3 11 



2Q 



Q 2 



+ 



ln (l-x+)(l-x+)g | (HQ 2 /Q 2 T/ 



Q 2 



+ 



(Qt)- 

1 



yr+ef x l x l 

1 dx 



+ 



+ S(Q 2 T ) 



71 

- 4 + T 



T^e-v x 2 x 2 

5H(x° l , x^ 1 



5H(xi , 
&ff"(:r* , x 2 ) 



1 >»£ 



T 



+• 



- 8H{x\ , X th 



xix* 2 

T 

x\x 2 



SH{x\ , xl) 



2C r 



x 



- A T P qq (z) 5H , x{] + f - A T P qq (z) 5H (x\ , ^ 

Jx\ Z \ Z J J x o Z V Z 



(4-23) 



To completely isolate the growth of the form ~ 1/<5t as Qt decreases from the third and 
fourth lines of (4-23), one more step is required. We rewrite the integral in the third line as 



dxi 



T^eV x l x \ 

1 



+ 
1 

dx\ 



5H(x 1: x* 2 ) --SH(x° 1: x° 2 ) 



X\X 2 



y^ev x l x l 



5H( Xl ,x* 2 ) - 8H{x\ , x%) 



dxi 



X\ — x\ 



-0 

„0\ 1 



6H( Xl ,x»)-±-6H(xl,x») 
x i 



I 



+ 



dxi 5Hi + 



dx\ 



X\ — x{ 



5H{x x , x l) — - 5H{xl , xl) 
x i 



2C f 



1 dz 



z 



A T P qq {z) 6H[^,xl)- 8H{x\ , xl) ( - A + In 



x\ 



(4-24) 
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introducing the following shorthand notation for the integral that vanishes for Qt = 0: 

1 dxi 



I dx 1 5H 1 = I 



^r e y x i x i 



5H(x\ , x* 2 ) 



X\X 2 



8H{xl , x° 2 ) 



dxi 



X\ — x\ 



X 







5H( Xl ,x u 2 )^- dHtxl , x° 2 ) 

X\ 



We also rewrite the integral in the fourth line in the same way, using 

dxo 



J dx 2 5H 2 = J 

-I 



1 



5H(x* 1 ,x 2 )- 5H(x\,xl) 



x\x 2 



dx-) 



x 2 - x u 2 



8H{x\ , x 2 ) ^ - SH{x\ , x° 2 ) 
x 2 







dxiSHi 



y—-y 



In this way, (4-23) finally becomes 

A T da R ( S )+v+CT 2?r „ 1 
" ^M,,, =7yi G (°) cos ^ ~ & ~ **) ~ 



dQ 2 dQ^dyd(p kl Q 
2 



6 



1 + 



x 



Qt 



chu 6Hi + / dx 2 SH 2 ) + 2 5H{x\ , x° 2 ) In 
J J Qt 



2 Qj 1 
(1 - x°)(l - x°) 



+ 8H{x\ , x° 2 
1 



+ 



( HQVQt 

Q 2 



+ 



(<£)- 



(-8 + 7t 2 )5(Q 



+ (J(Q^) In -\ 



x 



x- 
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(4-26) 



(4-27) 



Setting e = 0, the contribution to (4-8) from the nonsingular part of A T 1Z, (3-8), is easily 
calculated as 



A T da(x 1 P l , x 2 P 2 ) R{NS) 2tt h 



1 



dQ 2 dQ^dyd<p kl 



= -7^ G (°) cos ( 2 ^i - & - fa) 7yT ln 1 + ^7 



g 



g 



T 



Q\ 



Q 2 



X 



5(xiX 2 — X\X 2 — x 2 x\ + r). 



Thus, the contribution to the hadronic cross section, Z^g^r in (4-6), is simply 
A T da-n^NS) 



dQ 2 dQ\,dyd(j)k 1 



^ G H (0) cos (20 fcl - <p x - 2 ) ^ ^ (l + ||) 
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^/r+ey El 30 1 X\X 2 J ^/r+e~y -^2 %2 X^X2 



(4-28) 



Adding (4-27), (4-28), and the tree-level contribution in 4 dimensions [see (4-3)-(4-6) and 
compare with (4-21)], 



and using 



2tt 



Q 2 G//(0) 167r 2 iV c 5 V Q 2 J vr 



47ra\ a 



c F , 



we reach the final expression for the differential cross section (4-1) for the tDY (2-1) as (4-6), 
including the 0(a s ) QCD mechanism in the MS factorization scheme: 



A T da 
dQ 2 dQ 2 T dyd(j) kl 

Here, we have 



cos(20 fcl -0 1 -0 2 ) [A T X(Q 2 T ,Q 2 ,y) + A T Y(Q 2 T ,Q 2 ,y)}. (4-29) 



jV = 



4ira\ 1 



a 



16n 2 N c S \Q 2 J 3 3iV c ,SQ 2 ' 

AtX contains all terms that are singular as Q\ — > 0, behaving as Q^ 2 x (hKQVQr) or 1) 
or 

5(Qr)» and A T Y is the remaining "finite" part (containing less singular terms). Writing 
A T X as 

A T X (Q 2 T ,Q 2 ,y) = A T X® (Q 2 T , Q 2 , y ) + z\ T X« (Q 2 , Q 2 , y), (4-30) 

we have 



ArX<°>(Q* , Q 2 , y) = , x° 2 , 5(Q 2 T ), 

A T X^(Q 2 T ,Q 2 ,y) 

3 



(4-31) 



+ 



X 



2tt 
1 



/ln^QVQi; 

V Q 2 T 



+ 



(Q 2 t)- 



+ (tt 2 - 8) 5{Q 



5H(x»,x u 2 ;pi F ) 



my 



+ 5(Q 2 T )\n^ 



£^A T P qq (z)5H(^,x 2 ; ^+J^A T P qq (z)5H^xl^; ^) J, 



(4-32) 
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and 




x\x 2 

f 1 dx 2 




5H{x\ , x 2 ;n 2 F ) — 



x\x 2 



T 



Qt 



1 



/ 



dx\ 8 H\ + 



dx2 SH 2 + 5H(xi , x 2 ; up) In 



(l-x»)(l- a O)_ 



(4-33) 



where we have recovered the notation of (2-3) for the factorization scale with the replacement 
/i 2 — > n F , so that /jl f in (4-32) and (4-33) denotes the MS factorization scale. We have also 
made explicit the dependence of the QCD coupling constant a s (/i 2 ? ) on the corresponding 
MS renormalization scale /ir. When both of the colliding hadrons in (2-1) are polarized 
along the x-axis, we have 0i = 4> 2 = [see (2-12)]. Note that this result is invariant under 
the replacement y — > —y, as it should be. 

Equation (4-29) with (4-32) and (4-33) was first obtained in Ref. 23), but its derivation 
was not described in detail there. Here we give a complete derivation, explaining all the 
necessary techniques to deal with the complications in the D-dimensional calculation involv- 
ing transverse degrees of freedom. We note that, integrating (4-29) over Q\, we obtain the 
cross section A T da / dQ 2 dyd(j) kl in the MS scheme, and the result coincides with that found 
in Ref.20). In Appendix B, we also report the expression for A T dcr/dQ 2 dQ 2 ^d(f)k 1 , integrating 
(4-29) over y. When we integrate (4-29) over both Q\ and y, the result coincides with the 
corresponding total cross section associated with the partonic cross section (3-27), which was 
calculated in the last section. 

For Q\ > 0, Z\ T X (0) in (4-32) vanishes. Also, the terms proportional to 5{Q\) in (4-32), 
including those associated with the + distribution, ([\n n (Q 2 /Q%)\/Q\) + , do not contribute. 
The cross section (4-29) in this case is of 0{a s ), and it gives the LO QCD prediction for the 
spin-dependent cross section of the tDY in the large Qt region: 



A T da LO 



dQ 2 dQ T dyd(j) kl 



jV cos(20 fcl - 0! - 2 ) [a t X^ (Q2 , Q\ + A T Y(Q 2 T , Q 2 , y)] 



(4-34) 
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§5. Transverse-momentum resummation: general formulation 



The (^-differential cross section of the tDY including one- loop corrections, (4-29), con- 
tains the singular part A T X of (4-30)-(4-32) that grows as ~ a s \n(Q 2 / Q T ) / Q 2 T and ~ ct s /Q T 
as Q\ — > 0. The calculations in the last section indicate that these large contributions as- 
sociated with the singularities at Q 2 T = are induced as the recoil effects from the emission 
of a soft and/or collinear gluon. Actually, corrections of this type appear at each order of 
the perturbative calculation and become very large for Q T <C Q. Also, in the higher-loop 
corrections, the radiation of the soft and/or collinear partons produces the 1/e poles due to 
the IR divergences. After the factorization of the collinear singularities is accomplished in 
terms of the parton distribution functions, the remaining IR divergences are found to cancel 
when all the diagrams at the same order in a s are combined. The Kinoshita-Lee-Nauenberg 
(KLN) theorem 15 ) ensures that the corresponding cancellation of the IR divergences occurs 
between the real-emission and virtual-loop contributions, due to the relevant soft-collinear 
radiations. However, when the final dilepton of (2-1) is kinematically constrained to have a 
small transverse-momentum Qt, the real emission of the accompanying radiation is strongly 
inhibited. In this case, the KLN cancellation of the corresponding IR divergences still occurs, 
but it is "incomplete", because the phase space for the multiple real-emission in the final 
state is strongly restricted for small Q T ; specifically, all 1/e poles still cancel (at the leading 
power in l/Q) but the remainder depends in a singular manner on the small parameter 
Qt that constrains the phase space. These remnants of the incomplete cancellation at the 
boundary of phase space actually produce a series of enhanced logarithmic contributions, 

a r M m iQVQ 2 T) ( m = 2 n-l,2n-2,...,l,0) (5-1) 

at n-th order in a s . When only soft gluons are radiated and the sum of their transverse 
momenta balances Qt, the contributions corresponding to m = 2n — 1, as well as the 
contributions with m = 2n — 2, 2n — 3, . . ., arise as the IR-finite piece. When collinear gluons 
are also radiated, contributions with m = 2n — 2, 2n — 3, . . . arise.*) These so-called "recoil 
logarithms" make the fixed-order perturbation theory invalid for Qt <C Q, and have to be 
resummed to all orders in a s in order to make a reliable prediction of the cross section at 
small Qt- 

The corresponding resummation, the "QT-resummation," can be treated on the basis of 
the general formulation of Collins, Soper and Sterman (CSS). 32 - 1 After reviewing the CSS 
formalism, emphasizing its universal structure, we use 23 -* it to perform all-order resummation 

*) The higher-order contributions in (5-1) can be produced by the radiation of semi-hard gluons which 
have Qt as net transverse momentum. 
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of the recoil logarithms for the tDY cross section, up to NLL accuracy, which corresponds 
to completely summing the first three towers of logarithms, i.e., a™ln m (Q 2 /Q T )/Q T "with 
m — 2n — 1, 2n — 2, and In — 3, for all n. We also discuss various kinds of elaborations of our 
NLL resummation formula beyond the original CSS form, on the basis of recent developments 
in the QT-resummation formalism. 34 )' 35 )' 44 )' 45 ) Combining the resulting NLL-resummed cross 
section with the leading-order (LO) cross section (4-34) in a consistent matching procedure, 
we obtain the "NLL+LO" cross section of the tDY, which has a uniform accuracy over the 
entire range of Qt- 

5.1. Collins- Soper-Sterman (CSS) resummation formalism 

First, we explain the Qr-resummation formalism in a general form, such that it is not 
restricted to the present case of the tDY (2-1). For this purpose, we consider the process 

h^P,) + h 2 (P 2 ) - F(Q 2 , Q 2 T , y; 0) + X, (5-2) 

where the collisions of the hadrons hi and h 2 with momenta Pi and P 2 produce a system of 
non-strongly interacting final-state particles, F, carrying total momentum Q, total transverse 
momentum Qt and rapidity y. The additional variable <fi denotes the possible dependence 
on the kinematics of the final-state particles in F, such as individual transverse momenta. 
For the case of the tDY azimuthal angular distribution discussed in §4, we have <fi — > 0^, 
the azimuthal angle of the final lepton. In general, hi i2 may be polarized hadrons, but here 
their spins are suppressed for simplicity. 

The CSS resummation formalism yields the cross section for this process, which is appli- 
cable to the small Qt region as well as the large Qt region. The corresponding CSS formula 
possesses the following general structure, decomposed into two types of terms: 

da^ = da^ da^ 
dQ 2 dQ 2 T dyd<p ~ dQ 2 dQ 2 T dyd(/) + dQ 2 dQ 2 T dyd<p ' ^ ' ' 

Here, both terms on the RHS are expressed as convolutions of the corresponding partonic 
cross sections and the parton distributions, which are formally similar to the RHS of (2-3), 
but the partonic cross sections participating in (5-3) are more sophisticated than those 
calculated in §4. The first term in (5-3), do^^ j dQ 2 dQ 2 T dyd§, can be evaluated by resum- 
ming the "singular" terms, like AtX of (4-30)-(4-32), to all orders in a s . In particular, 
da( Tes:> /dQ 2 dQ' T dyd(f) collects the logarithmic contributions (5-1), which dominate the cross 
section at small Q T - The second term, da^ 11 ^ / dQ 2 dQ' T dyd4> ) is the "finite" component that 
is not associated with the logarithmically-enhanced contributions (5-1), and thus it can be 
computed with fixed-order perturbation theory. Formally, we can define the corresponding 
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finite component analogously to AjY in (4-29) and (4-33); i.e., da^ 11 ^ /dQ 2 dQ^dyd(J) is ob- 
tained as the contributions less singular than \jQ\ or 8{Q 2 r) as Q\ — > among the terms 
in the cross section do j dQ 2 dQ 2 ^,dyd<f), which is calculated up to an appropriate order in a s . 
In (5-3), the first term dominates the cross section d(r( css ) /dQ 2 dQ 2 T dyd<p for Q T < Q, and 
thus the second term is negligible in this region. But the second term becomes important 
when Q T ~ Q. 

In the CSS resummation formalism, the resummation of the logarithmic contributions 
(5-1) is carried out with the "6-space" resummation approach, introducing a 2-dimensional 
impact parameter b that is the Fourier conjugate of the transverse momentum Q T . As noted 
above, multiple-gluon emission induces (5-1) as its IR-finite piece, and the corresponding con- 
tributions can be resummed most straightforwardly in the b space, to all orders in the per- 
turbation theory. 26 - 1 After the resummation, the 6-space cross section is Fourier-transformed 
back to the Qt space. In this way, the first term in (5-3) is obtained as 32 ) 

*x (res - } i r „, „ „ r „ b 



dQ 2 dQ 2 r,dyd(f) 



i f r°° h 

— Jd 2 be^W(b;Q,y, ( p) = J db - J (bQ T )W(b; Q,y, 0), (5-4) 



where Jo(bQ T ) is a Bessel function, and the 6-space representation W (b; Q , y , 0) has the 
general structure [using the notation of Ref. 34) and the hadronic variables (4-7) also for the 
present case (5-2)] 

W(b; Q, v,<f>)=Y\ f dx x f dx 2 da(L ^ Q \ {Q 2 - Xl x 2 S)5 (y - \ In ^ 
^io io # V 2 x 2 



x{Cji®fi /hl ) (x lt ^ (C ]k (g) f k/h2 ) (x 2 



IP 



^ d ^^W^Q,4,4), (5-5) 

j 



with 



Wj(b; Q, Xl ,x 2 ) = e s ^ Q \C 3t ® f i/hl ) (x lt |) (C ]k ® f k/ha ) (x 2 , |) . (5-6) 

Here, the subscripts j, i and k can be q, q or g, including the flavor degrees of freedom, 
cr^ )(Q 2 ) is the lowest-order cross section (integrated over Q 2 ^)*^ for the parton-level process 
j + j — > F(Q 2 ), and b = 2e~" /E , with ^ E being the Euler constant, has a kinematical 
origin. The quantity Cji is the coefficient function, and fi/h(x,fi 2 ) is the unpolarized or 

*) For the spin-dependent cross section in tDY, j = q is relevant, where c'fo^Q 2 ) is the LO cross section 
for the qq annihilation process q + q — > ll(Q 2 ), which is equal to Atct™ calculated in §3 [see (3-25), (3-27)]. 
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polarized distribution function for a parton i inside the hadron h. The symbol ® denotes 
their convolution, defined by 

r 1 d z 

(Cji® fi/h) (x,fi 2 ) = —Cjiiz^asifi^fi/hix/z,^ 2 ). (5-7) 

J X Z 

The large logarithmic corrections are resummed into the Sudakov form factor e Sj ( b '^ in 
terms of the exponent 

Sj(b, Q) = - A 3 {a s {n 2 )) In ^ + B^k 2 )) , (5-8) 

with the perturbatively calculable functions Aj(a s ) and Bj(a s ). Here, for Q ^> 1/6, the 
leading contribution from the Aj(a s ) term is enhanced in comparison with that from the 
Bj(a s ) term due to an explicit logarithm ln(Q 2 / k 2 ) in the integrand. Other large logarithms 
are also implicit in the integration over k 2 for both terms. In the form factor e Sj( - b '®\ the 
Aj(a s ) term of (5-8) represents exponentiation of large logarithms due to soft radiation, 
while the Bj(a s ) term represents exponentiation of large logarithms due to flavor-conserving 
collinear radiations. 44 -* 

The quantity Wj(b; Q, x±,x 2 ) in (5-6) depends on Q 2 only through the Sudakov exponent 
(5-8). This simple structure is a consequence of the fact that Wj(b; Q,Xi,x 2 ) is obtained as 
the solution of the following evolution equation: 32 ) 



d 

Wj(b; Q,xi,x 2 ) 



dlnQ 



^Ma s (K 2 )) + BjMQ 2 )) 
bl/V K 2 



W j (b;Q,x 1 ,x 2 ). (5-9) 



Moreover, the large logarithm of Q 2 b 2 , arising in the integral J b 2^ b2 (dK 2 / K 2 )Aj(a s (K 2 )) of 
(5-9) when b ~ l/<5r ^ is controlled by the renormalization group (RG) equations 

governed by certain anomalous dimensions. For j = q, which is relevant to the DY and vector 
boson production, the corresponding RG equations, as well as (5-9), have been derived on the 
basis of the factorization property of the quark form factor in QCD into "hard," "soft," and 
"jet" factors to all orders in a s , and also, exploiting the renormalization property of those 
individual factors, defined as matrix elements of gauge-invariant operators (see Refs. 32) and 
48)).*) 

The solution of (5-9) is generally expressed as 

W J (b;Q,x l ,x 2 )=e s ^ W J (b;b /b,x 1 ,x 2 ), (5-10) 



*) In this approach, A q and B q appearing on the RHS of (5-9) are related to another set of two evolution 
kernels, denoted as "if" and "G" in Ref. 32). K and G, respectively, obey the RG equations governed by the 
so-called cusp anomalous dimension. This relation allows us to confirm that A q and B q contain no h\Q 2 b 2 
contributions and indeed take the forms of power series in a s . 
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using the Sudakov exponent Sj(b,Q) of (5-8). With the boundary value Wj(b;bo/b,x\,X2) 
specified, (5-10) determines the complete behavior of Wj(b; Q, xi, x 2 )- For this purpose, 
we note that when Qt <C Q, and for the relevant region b ~ \/Qt, large logarithms of 
Q 2 b 2 ~ Q 2 /Qt indeed arise in the Sudakov factor of (5-10) via the integral of (5-8). On the 
other hand, the boundary value Wj(b; bo/b, xi, X2) depends explicitly on only one distance 
scale, namely b, and thus, for b <C 1//1qcd, its b dependence is calculable by the customary 
perturbation theory. Specifically, Wj(b; bo/b, x\, X2) is determined by the singular component 
of the fixed-order differential cross section, which has a form analogous to that of AtX 
in (4-30)-(4-32) obtained in §4. Setting Q = b /b and the factorization scale as up = 
bo/b in formulae like (4-30)-(4-32), the Fourier transformation of the corresponding singular 
component from the Qt space to the b space yields a result with the structure 

Wj(b; b /b,x 1 ,x 2 )=J2(C ji ®f i/hl ) (C ]k ® f k/ha ) (x 2 ,^j , (5-11) 

which is factorized into two parts, corresponding to the incoming hadrons, hi and h 2 . This 
reflects the collinear nature of radiative corrections associated with the singular terms in 
the relevant region, Qt ~ 1/6; i.e., the cross section receives logarithms of Mqcd only from 
collinear radiation, which can be treated with the DGLAP evolution of the parton distribu- 
tions fi/h(x,bl/b 2 ) associated with each hadron, and Cy t of (5-11) represents the remaining 
perturbative corrections as a power series in a s . Therefore, combining the evolution equation 
(5-9) with (5-11) obtained from fixed-order perturbation theory, the general formula (5-6) 
follows. 32 ^ This analysis also demonstrates that (5-6) is accurate in the region b <C 1//1qcd- 
We express the expansions of the functions Aj(ct 8 ) and Bj(a s ) in (5-8) and the coefficient 
function Cji(z, a s ) in (5-6) in powers of a s as 

00 

^) = E(S)"4"'. 

n=l 

00 

b,<-) = E (£)"*?". 

n=l 

00 

C 3l (z,a s ) = 6*6(1 -z) + J2 (^f^) • (5- 12 ) 

n=l 

The coefficient yields the LL resummation, {A^ , , C^} give the NLL terms, 
{A^ , B^ , C^} give the NNLL contributions, and so on. In particular, among the coeffi- 
cients {A^\ A {2 \B {1 \C W } that are necessary for the NLL resummation, , A (2) , B^} 
are known to be independent of the final states (see Refs. 34), 44) and 35)). For the DY 
process, (5-2) with F — 11, which we consider in this paper, and also for other processes, 
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such as W and Z boson production, these coefficients with j — q are necessary, and they are 
given by 



A^=2C F , Af = 2C F K, B^ = -3C F , (5-13) 



wr 



ith 29 ) 



where Cq = N c , and Nf is the number of QCD massless flavors. The result (5-13) can be 
derived directly by evaluating certain loop diagrams which represent the evolution kernel of 
( 5 . 9 )_32),48) Apparently, for A { q 2 \ this requires a two-loop calculation. Nevertheless, its value 
is known to be independent of the process as well as the spin. In fact, we can confirm that 
this is indeed the case by using a relation 29 ^ between and the usual two-loop DGLAP 
kernels for the parton distribution functions. *•* The parton distributions in (5-6) obey the 
DGLAP equation, 14 ) 

f^fi/K^f) = (? as(/i2) ) M*'^ ( 5 - 15 ) 



j 

with the corresponding DGLAP kernel &ij(z, a s (fi 2 )) given by a power series in a s : 

^(z, a s ) = ^\z) + g 0>\)\z) + ■■■. (5-16) 

The behavior of &ij(z, a s (fi 2 )) for z — > 1 is dominated by soft gluon emissions, and it 
is diagonal in Indeed, the coefficients A^ and A^ are, respectively, related to the 
dominant large- z behavior oc 1/(1 — z) + in the one- and two-loop terms of (5-16) for emission 
of gluons from a quark: 

The large z behavior (5-17) is universal to all DGLAP kernels associated with the twist-2 
quark distributions, i.e., to the DGLAP kernels for density distribution q(x,fi 2 ), helicity 
distribution Aq(x, fi 2 ), and transversity distribution 5q(x,fi 2 ). We also note that the term 
oc 5(1 — z) in the one-loop kernel &qq(z) is also universal, and its coefficient determines the 



value of Bg in (5-13). 34 ** The similar large- ,2 behavior 

1 + <?V)( Z \ ^ 2Cg + ^1 2CgK 
99 y~> ^ 2vr 99 y } (1 - z)+ 2tt (1 - z)+ 



*) The coefficients and A q 2 ^ are given by the one- and two-loop terms of the universal cusp anoma- 
lous dimension, respectively. The cusp anomalous dimension plays a role 73 ) in the evolution of the quark 
distribution functions for large x. 
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and the universal coefficient of the 5(1 — z) term in the one-loop kernel ^gg(z) lead to 49 - 1 

4 1) = 2C G , Af=2C G K, 5« = -47rA), 
with K of (5-14). Here, f3 is the first coefficient of the QCD (3 function, and is given by 

A) ^ • (5 18) 

Recently, the three-loop term ^^\z) in (5-16) for the unpolarized parton distributions has 
been calculated, and A^p has been extracted from its z — > 1 behavior. 50 ) We note that 
(n > 1) and, thus, the entire function Aj(a s ) are actually process independent. As for the 
other coefficients in (5-12), (n > 2) and Cj^ (n > 1) depend on the process. 34 ) General 
expressions for the coefficients and Cj]\ including the process-dependent pieces, are 
derived in Ref. 34) in the MS scheme.*) It is possible to transform the CSS resummation 
formula given in (5-4)-(5-6) into a new form, so that the process independence of the building 
block of the resummation formula is maximal (see Ref. 44)). 

The logarithms of Q 2 b 2 contained in the Sudakov factor e s ( b ' Q \ which are large in the 
region b ^> 1/Q, can be made explicit. Here we have suppressed the subscript j of the 
Sudakov exponent for simplicity. These contributions can be organized within a systematic 
large logarithmic expansion according to consistent order counting, where a s \ia(Q 2 b 2 ) is 
formally considered of order unity. 35 )' 45 ) Substituting the expansions (5-12) and the explicit 
form of the running coupling constant, 

Gs[K } A)M«V^cd) $1bV/4cd) + ' ' ' ' l "' i9j 
with p of (5-18) and 

_ 17C 2 G -5C G N f -3C F N f 

Pl ~ 24^ ' (5 ' 20) 

we perform the k 2 integral in (5-8) and organize the resulting exponent S(b,Q) in terms of 



Q 2 b 2 

A = /3 a: s (/4) In -T2- = (3 a s (n 2 R )L, (5-21) 



where L plays the role of the large logarithmic expansion parameter in the b space. Consid- 
ering as A ~ O(l), S(b,Q) can be systematically expanded as 

-. oo s , 2 \ \ n—1 

S(b, Q) = -i^ M°)(A) + ft«(A) + J2 fcW(A), (5-22) 

*) (n = 1,2, . . .) and B^ are independent of the factorization scheme, but Bq 1 ^ (n > 2) and c\™\z) 
(n > 1) depend on the factorization scheme (see e.g. Ref. 35)). 
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where h^ n \\) (n = 0, 1, . . .) are functions of A, which vanish at A = 0, as h^ n \0) = 0. Then, 
it is straightforward to find 

AW 



hW(\) 



2n(3 2 



[A + ln(l-A)], 



(5-23) 



A) 



BW 



ln(l - A) 



Air 2 ? 2 



Q 2 



A< 2 >- 2tt / M (1) ln^ 
/4. 



A 



1- A 



+ ln(l-A) 



(5-24) 



and /i^ 2 ^(A) has a similar form, which can be found in, e.g., Ref. 35); ft/ 2 ) (A) involves 
.B^ 2 ) and the third coefficient of the QCD (3 function, fo, in addition to the coefficients 
appearing in (5-24). From (5-22), one can see that the Sudakov factor e s ^ b '^ is generically 
the exponentiation of the logarithmic terms a™L m with m < n + 1 at each order n: The first 
term, h^/a s , in (5-22) collects the LL terms, a n s L n+1 (n > 1), and the second term, /^^(A), 
collects the NLL terms, ct n s L n (n > 1). Also, the third term, (a s /2%)h( 2 \ controls the NNLL 
contributions, (n > 1), and so forth. Note that the ratio of two successive terms in 

(5-22) is formally of 0(a s ), and (5-22) is organized as a systematic power series in the small 
expansion parameter a s (/i 2 R ), similarly to the customary perturbative expansions. 

5.2. Matching with the fixed-order a s calculation 

To bring the CSS resummation formula (5-4)-(5-6) into contact with the fixed-order 
calculation of the cross section, we formally expand W{b; Q, y, 0) in (5-4) in terms of a s (fi 2 R ), 
and perform, order-by-order, the Fourier transformation to Qt space. Here we consider 
the expansion of if (b; Q , y , 0) up to 0(a s ). Using (5-22), we can immediately expand the 
Sudakov factor e 5 ^'^ in terms of a s (fi R ) as 



S(b,Q) 



1 - 



Ml (\ A W L 2 + bWl]+ 0(a 2 ), 
2n \2 J 



(5-25) 



with L = \n(Q 2 b 2 /ft 2 ,)- Also, the b dependences of the coefficient functions Cji and the parton 
distributions fi/h in (5-5) can be expanded perturbatively in terms of a s (fi 2 R ); specifically, 
(5-12) and (5-15) yield 



2tt 



{z,a s (bl/b 2 )) = 6^6(1 -z) + ^^cW(z) + 0(a 2 s ), 
fi/h(x,b 2 /b 2 ) = f i/h (x,fi 2 F ) 



2tt 



(5-26) 
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Substituting these expansions into (5-4), we encounter an integration of the type 

/„ = / (W-MbO-r) I hi— ^- I . (5-27) 



o 



b 2 



This can be calculated by taking the i] — > limit of the similar Fourier transformation of 
(b 2 Q 2 /bffl [see (6-14) and (C-3) below]. For the present case, we need the results 

' In Q 2 /Q 2 T ^ 



1q = 8(Qt) > ^i — 



1, = -2 



(5-28) 



where the + distributions are defined as (4-19), and we get 



do (xcs 



dQ 2 dQ 2 T dyd(j) S 



dall Q) (Q 2 ) 



2tt 



HQVQl) 



KQt) fj/h! (xl,fJ? F ) fj/h 2 (xl,fi 2 F ) 



+ B 



(i) 



Qt J+ 1 (Q¥)+ 

Dql?®/^) 04/4) 



fj/h f]/ h2 (A,n 2 F ) 



+ 



E^i^^JW^D/j/^ 04/4) 



(5-29) 



Here, we identify the structure characteristic of the singular component of the fixed-order 
cross section, like that observed in A T X of (4-30)-(4-32). In particular, the matching of (5-29) 
with the corresponding fixed-order cross section completely determines the coefficients 

and and the LO cross section <Xjx )(Q 2 ) for the parton-level process j+j — > F(Q 2 ); 
for example, for the case of the spin-dependent cross section of tDY, the matching with 
(4-30)-(4-32) gives results identical to (5-13) for and and also yields 



C$(z) = C$(z) =Cf(j- 4) 5(1 - z), C$>(z) = Cg\z) = 



(5-30) 



in the MS scheme and 

ld^[ 0) (Q 2 ) _ldafl Q) (Q 2 ) 



= (? cos (20 fcl - 01 - 2 ), 



S d<p S d<p 

while all other quantities associated with the gluon index j — g vanish. 



(5-31) 
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The expansion (5-29) can be extended to include higher orders in a s . The corresponding 
expansion up to 0(a 2 ) was carried out for Qt > in unpolarized DY, 30 ) and the result- 
ing 0(a 2 ) term involves the coefficients and B^.*^ The matching of that result with 
the corresponding fixed-order cross section confirms the universal value of in (5-13) 
and also determines for unpolarized DY [see also Ref. 34) for results in other unpo- 
larized processes of the type (5-2)]. This fact also implies that when we use the resummed 
component (5-4) up to NLL accuracy, da^ 11 ^ j dQ 2 dQ\dyd<\) in (5-3) should be taken as the 
finite component of the fixed-order a s cross section, and the sum of those two components 
gives the cross section da^ css ^ / dQ 2 dQ 2 ^dyd(j), which is exact up to 0(a s ) when expanded 
in powers of a s . Similarly, when we use the resummed component up to NNLL accuracy, 
dcr( fin ') / ' dQ 2 dQ\dyd(j) should be taken as the finite component of the fixed-order cross section 
to 0(a 2 ), so that we obtain da^ css ^ /dQ 2 dQ' T dyd(f), which is exact up to 0(a 2 ). Generalizing 
these considerations, the second term, da^ 11 ^ /dQ 2 dQ 2 ^dyd(j), in (5-3) should be determined 
by 

da^ _ da 
dQ 2 dQ 2 ^dyd(j) dQ 2 dQ' T dyd(f) 

where the notation (• • • )| fo represents the expansion of the quantity (■ • • ) in powers of 
up to a given fixed order; i.e., the first term on the RHS, da /dQ 2 dQ^dyd(f)\ {o , is the cross 
section that is computed by truncating the customary perturbative expansion at a given 
fixed order in a s (fi 2 R ). As noted above, when we consider the resummation at the NLL 
level, (•••)!&> is the expansion up to 0(a s (fj, 2 R )). In this case, the second term of (5-32), 
da^ /dQ 2 dQ 2 T dyd(p\ io , is given by (5-29), where the coefficients A™ and 5« in (5-12) 
and ^(lo)^ 2 ) > determined similarly to (5-30) and (5-31), are substituted. When we 
consider the resummation at the NNLL level, (• • • )|fo is the expansion up to 0((a s (fi 2 R ) 2 ). 

The matching procedure for the finite component, represented by (5-32), guarantees that 
cfo-( css ) / 1 dQ 2 dQ\,dyd4> in (5-3) retains the complete information of the perturbative calcula- 
tion up to the specified fixed order and, at the same time, incorporates the resummation 
of the logarithmically enhanced contributions to all orders. The fixed-order truncation of 
da^ css ydQ 2 dQ^dyd(j) exactly reproduces the fixed-order cross section obtained in the cus- 
tomary QCD perturbation theory. In this sense, using the matching of (5-32), the resummed 
and fixed-order components are consistently combined in da^- css ^ / 1 dQ 2 dQ\dyd(j) of (5-3) with- 
out any double counting. 

*) The 0{a 2 s ) term involving is proportional to 8(Q T ), similarly to the 0(a s ) term involving 
in (5-29), and it vanishes for Qt > 0. 



fo 



da^ res -> 
dQ 2 dQ 2 ^dyd(f) 



(5-32) 



fo 
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§6. Systematizing the b-space resummation formula 



In the b integration in the resummed component (5-4), we can distinguish three re- 
gions: 26 )' 32 ) (i) < b < 1/Q, (ii) 1/Q < b < 1/v1 Q cd, and (iii) b > 1/A QCB . These are 
relevant to the kinematical regions Q T > Q, ^Iqcd <C Qt ^ Q, and Q T < /Iqcd, respec- 
tively. We now go into further detail concerning the behavior and the physical content of 
the resummation formula in each of these three regions. 

6.1. Region (i): < b < 1/Q 

In the region (i), < b < 1/Q, the resummation is irrelevant and can be truncated at 
low order in a s . Thus here, the customary fixed-order perturbation theory can be used. The 
resummed component associated with this region should be given by (5-29) with Qt > Q to 
good accuracy. Apparently, in general, if the functions A, B and C of (5-12) and the QCD 
(3 function are evaluated up to order a™, then the corresponding resummed component is 
correct to order a™. 

6.2. Region (ii): 1/Q < & < 1/A QCD 

In the region (ii), 1/Q < b <C 1/Aqqb, all-order resummed perturbation theory can be 
used [see the discussion below (5-10)]. In this region, L in (5-21) can be large, but still we 
have, taking jiR ~ Q as usual, 

< A < A)a s (/4) In -S- ~ 1, (6-1) 

yl QCD 

so that the Sudakov factor e s( - b '^ with (5-22)-(5-24) can be expanded in a power series in A. 
In principle, the b dependence of the coefficient functions Cji and the parton distributions 
fi/h in (5-5) can be expanded similarly in powers of A (see §7). 

It is instructive to apply the corresponding expansion to (5-4) and (5-5) at LL accuracy, 
which means that we retain only the first term with in (5-22) for the Sudakov factor 
e s(b,Q) ^ anc j a j go t erms at the same level in the other factors. The coefficient function 
in (5-12), as well as the scale dependence of the parton distributions fi/h, appears at the same 
logarithmic level as i.e., at the NLL level, as seen in (5-29) (also compare (5-24) and 
(7-17) below). Thus, at the LL level, those effects can be ignored as (Cji <S> fi/h)(x, b^/b 2 ) — > 
fj/h(x,fip) in (5-5). As a result, the b dependence of W (b; Q , y , 0) in the integrand of (5-4) 
remains only in the Sudakov factor at the LL level [see (5-23)], 



e s(b,Q) = e M°)(A)/a s 




n=0 \ r{i s m=2 
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n 



n=0 



(6-2) 



where a s = a s (/i%). In the second line, (5-21) has been substituted, and the series are 
rearranged according to their powers in the large logarithm L for each order in a s , with 
von' 1 = (— / 4ir) n , and other expansion coefficients, Wn\ Wn \ ■ ■ •, which can be expressed 
similarly in terms of A^ and (3 . Using the expansion (6-2), the Fourier transformation 
relevant to (5-4) at the LL level is given by 



poo 1 

/ db b -J (bQ T )e s ^ = T^(Q 2 T ,Q 2 
Jo * 



?1=0 



(6-3) 



and it can be completely expressed in terms of X m (m = 2n, 2n— 1, . . .) in (5-27). Generalizing 
(5-28) for an arbitrary integer n (> 0), we get [see (C-3)], 



T = 

- L -"n 



drj 

qJ{Qt) ~ n 



5{Q 2 T ) - J2 



V m+1 fln m (Q 2 /Q 2 T ) 



m=0 
.n-1 



ml 



Q 



+ 



exp 



- 2 E 



r=l 



Cgr+l) 2r+1 
2r + 1 ' 



?7=0 



ln^ 1 (Q 2 /Q 



Q 2 - 



+-C(3)n(n-l)(n-2)(n-3) 



ln^(QVg|) 



+ 



(6-4) 



'T 

where g n = (d n /dn n ) exp [— 2 X^li C(2?" + l)r? 2r+1 /(2r + 1)]| =0 , C( n ) is the Riemann zeta- 
function, and the ellipses denote the terms involving (In™ (Q 2 / Q\) / Qt) + with m < n — 6. 
Thus, for /Iqcd < <5t < <5, corresponding to the region (ii) l/Q < 6 < 1/jIqcd, (6-3) yields 

A^\ n \n 2n -\Q 2 /Q 2 T ) 



n=l 



t^q 2 t ,q 2 )=j:^ 



-In 



47T 



Q 



+ 



(6-5) 



where the ellipses denote the terms involving In™ (Q 2 / 'Q\) / 'Q\ with m = In — 2, In — 3, • • • . 
This result demonstrates that (5-4) indeed resums the towers of logarithms (5-1) to all orders 
in a s . Summing up the leading tower in (6-5) over all n, we get 



T-HQl,Q*)- AmaMQVQl) 



27T 



Q 



exp 



A 0) 



a.. 



An 



\n 2 {Q 2 /Q 2 T ) 



+ 



(6-6) 



Substituting A^) of (5-13), it is seen that this result reproduces the Sudakov quark form 
factor in the so-called "double leading logarithmic approximation (DLLA)." 24 ) 
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It is worth noting that, introducing the Fourier transform [see (5-28)], 



u(k T ) = 



cPb 
An 



-ibk,T 



AVas 
Arc 



2tt V 



(6-7) 



the contribution from the leading tower of logarithms in (6-3) can be expressed as 

oo 1 

T- 

=tt J] - y d 2 A; 1)T • • • d 2 k n , T 5^ I Q T + J2 



n=l 







J An 


An 



v{K,t) 



(6-8) 



n=l 



Here, z/(/c^) in (6-7) represents the probability distribution for the emission of one soft gluon 
with transverse momentum k? in the final state, as (da 1 '^ 011 / dk 2 ^) / a total = f(A/r), with <r total 
the total cross section. 26 ) Indeed, terms with the same structure as (6-7) appear in (4-32) and 
(5-29). Therefore, the n-th order term in (6-8) represents the probability for the emission of 
n soft gluons in the final state as the product of the probabilities for n independent gluons. 27 ) 
Thus, the exponentiation of the LL term in the Sudakov factor e s ( b '^ in fact sums up the 
multiple soft-gluon emission probabilities to all orders. 

Note that in each term of the expansion on the RHS of (5-29), transverse-momentum con- 
servation is ensured by the delta function, 5^ (Q T + YH=i ^i,t)i an d thus Q T of the dilepton 
in the final state is provided by the recoil from the radiation of soft gluons. This automatic 
and proper treatment of transverse-momentum conservation is one of the most important 
features of the impact parameter 6-space approach. 26 > In fact, transverse-momentum conser- 
vation is particularly important when we approach the small Qt region. For example, the 
contribution for n = 2 in (6-8) gives, using (6-4), 



n 
2! 



] J d 2 k hT d 2 k 2 , T 5^ (Q T + k ltT + k 2 , T ) ^ 



u(k 1>T ) v{k 2 , T ) 



n 



a, 



2n 



\n\Q 2 /Q 2 T ) 



Q 



+ 



+ 



4C(3) 
(Q 2 t) + 



(6-9) 



where the first term corresponds to the n = 2 term given in (6-5). Inspecting the integration 
over the transverse momenta, k^ T and fc 2i r, on the LHS of (6-9), we confirm that the first term 
on the RHS comes from the "strongly-ordered" phase space as 5 <y2 \Q T + k^ T )9(k\ T — fcf T ) + 
(1 <-> 2). Then, the second term is the non-strongly-ordered contribution corresponding to 
the emission of soft gluons whose transverse momenta, fc ljT and k 2 ,r Qt), add vectorially 
to balance Q T . Although the second term is formally subleading in (6-9), contributions of 
this type, called "kinematic logarithms," 27 )' 51 ) dominate the cross section in the asymptotic 
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limit, Q ^> Qt ^Iqcd, where the DLLA Sudakov factor of (6-6) leads to suppression of 
the (formally) leading contributions. Here, the probability of having a parton-antiparton 
annihilation into a lepton pair with no emission of gluons of transverse momenta greater 
than a fixed value decreases asymptotically faster than any power of Q 2 , and events at 
small Qt may be obtained asymptotically only by the emission of at least two gluons whose 
transverse momenta are not small and add to a small value of Qt- All the non-leading 
"kinematic logarithms" are correctly taken into account by imposing transverse-momentum 
conservation as in (6-8) and (6-9), and this is naturally realized in the 6-space resummation 
approach with (5-4).*) Also, apparently, transverse-momentum conservation with the 6-space 
resummation approach is crucial in the case of smaller Qt corresponding to the region (iii). 

Calculations similar to those in (6-2)-(6-9) can be performed including the NLL terms and 
higher logarithmic terms. It is straightforward to show that the NLL terms associated with 
of the Sudakov exponent (5-22) give rise to a series of contributions a™L 2n_1 , a™L 2ri ~ 2 , . . . 
to (6-3), and thus yield contributions a™ W n {Q 2 /Q T )/Q T with m = 2n — 2, 2n — 3, . . . to 
(6-5); similarly, the NNLL terms give rise to the contributions < \n m (Q 2 /Q 2 T )/Q 2 T with 
m = 2n — 4, 2n — 5, . . . to (6-5). Therefore, the resummation formula (5-4) at the LL level 
corresponds to summing up only the first tower of logarithms, a n s \\y n ~ x {Q 2 jQ\) / Q T exactly, 
while the resummation formula up to NLL accuracy allows us to fully sum up the first three 
towers, i.e., a n s \n m (Q 2 / Q 2 T ) / Q 2 T with m = 2n - 1, 2n - 2, and 2n - 3. 

The convergence of the resummed perturbation series for (5-4) in the Qt space is governed 
by a s \n 2 {Q 2 /Q't) rather than a s , and when Qt is small, a s \n 2 (Q 2 /Q' T ) will be large even 
for a s <C 1. As Qt approaches Q, we have ln(Q 2 /Q T ) <C 1, and the all-order resummed 
perturbation theory of (5-4) is smoothly extrapolated to the region (i). 

6.3. Region (iii): b > 1/^.qcd 

Next we turn to the region (iii), b > 1/tIqcd- In this region, A in (5-21) can be as large as 
1, and hence the Sudakov factor e s ^ b ' Q ^ cannot be expanded in terms of A, and the functional 
form exponentiating (5-22)-(5-24) has to be used as it is. Moreover, in this long-distance re- 
gion, nonperturbative effects become relevant, which could modify the resummation formula, 
(5-5) and (5-6). In fact, when b becomes extremely large, so that A — > 1, the resummation 
formula in the form of (5-5) and (5-6) breaks down. In this case, the functions (5-23) and 
(5-24) in the Sudakov exponent (5-22) are singular at A = 1, and this singular behavior is 
related to the presence of the Landau pole in the perturbative running coupling a s (/t 2 ) in 

*) There have been efforts to formulate a "(J^-space" resummation approach, which is basically organized 
according to the strong ordering of the phase space. Thus this approach requires a hard and sophisticated 
task of including effects of the nondeading "kinematic logarithms" (see Ref. 51)). 
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QCD. 

Therefore, to properly define the b integration in the resummation formula (5-4) for the 
corresponding long-distance region, relevant to Qt < ^qcd, a prescription to deal with this 
singularity is required. Such a prescription may be accompanied by the modification of 
(5-6) to complement the relevant nonperturbative effects. For this purpose, two types of 
approaches are commonly used in the literature. *** One is the so-called b* prescription, which 
was proposed by CSS. 32 ) In this approach, Wj(b; Q,x±,x 2 ) in (5-6) is replaced as 

Wj(b; Q,x u x 2 )^Wj(K; Q,x 1 ,x 2 ) F NP (b; Q,x 1 ,x 2 ), (6-10) 

where 

K= . b == , b lim ~0.5GeV-\ (6-11) 
V 1 + b / fe L 

Through this replacement, Wj(b; Q,Xi,x 2 ) in (5-6), which is accurate for b <C 1//1q CD , is 
smoothly extrapolated to the extremely large b region, while the singularity is avoided, since 
K < &iim in Wj(b*; Q : x 1: x 2 ). This corresponds to effectively taking into account the "freez- 
ing" of the running coupling in the long-distance region due to the nonperturbative effects. 
The function F NP represents possible nonperturbative effects, which may be interpreted as 
the contributions associated with the intrinsic transverse momentum of partons inside the 
colliding hadrons. F NP is normally taken as a smearing factor of Gaussian type in b. 32 ^ Thus, 
in the large b region of the integration over b in (5-4), Wj(b*] Q, xi, x 2 ) approaches a constant, 
Wjipwrn] Q,xi,x 2 ), and F NP acts as a damping factor. Many different parameterizations of 
F NP have been proposed (see Ref. 36) and references therein). 

Another approach 45 ) consists of an extension of the minimal prescription 53 ) proposed in 
the context of the so-called threshold resummation. This approach allows us to avoid the 
Landau singularity in a purely perturbative framework: Decomposing the Bessel function in 
(5-4) into two Hankel functions as 

J (bQ T ) = \[H^\bQ T ) + H$\bQ T )] , (6-12) 

we deform the 6-integration contour associated with H^\bQT) and H^ibQx) into the upper 
and lower half planes in the complex 6-space, respectively, and obtain the following two 
convergent integrals as |6| — > oo: 

J_(rcs.) r L 

+ j db b -H<»(bQTW(b;Q,y,4;). (643) 

*) See Ref. 52), for an attempt to incorporate the nonperturbative effects as power corrections to the 
RG equations that control the large logarithms arising in the integral on the RHS of (5-9). 
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The new contour C± is taken to be from to b c on the real axis, followed by the two branches, 
b = b c + e t with t G {0, 00} and < 9 < 7r/4. The constant b c here is arbitrary in the 
interval < b c < bi, where bi = (bo/Q)e 1 ^ 2l3 ° as ^R^ corresponds to the position of the 
singularity on the real axis as the solution of A = 1. Equation (6-13) provides us with 
a (formally) consistent definition of a finite 6-integral for the resummation formula within 
a perturbative framework with which, unlike in the above case of the ^-prescription, no 
extra cut-off parameter &i im is required in this contour deformation prescription. In this 
paper, we employ (6-13) to deal with the Landau singularity. (Also, in this case, possible 
nonperturbative effects in the large b region can be included, as seen in (9-3) below. *)) 

We note that the contour deformation used to obtain (6-13) can be performed "safely" 
order-by-order in a s (fi 2 R ), by expanding W(b; Q , y, 0) in the integrand in powers of a s . 
The above discussion concerning the regions (i) and (ii) implies that for (5-4) in the region 
Qt ^> ^Iqcd we can define W(b; Q , y,<p) in the integrand as a power series in a s (fi 2 R ), which 
can be organized similarly as the integrand of (6-3). In this case, (5-4) can be completely 
expressed by X n in (5-27), whose b integration is well-defined and yields distributions like 
(5-28) and (6-4). It is straightforward to show that we can perform the above contour 
deformation into C± with (6-12) for the b integration in X ra ; in fact, on the LHS of the 
relation 



the corresponding contour deformation can be performed using Cauchy's theorem; i.e., the 
contributions to the integral from the contour in the limits 6 — » and |6| — > 00 are found to 
vanish using an appropriate analytic continuation for 77. After such a contour deformation 
is performed for each term in the series of (5-4), the corresponding series can be summed 
up under the integrand of the deformed ^-integration. This yields W(b; Q , y,4>), with the 
Sudakov factor exponentiating (5-22). This final form coincides with (6-13) exactly, and is 
applicable to the region Q T < /Iqcd as well as the region Q T ^> Aq CD . Therefore, the 
choice of contours in (6-13) is equivalent to the original contour in (5-4), order- by-order in 
a s (/i R ), and also extends the applicability of the formula even into the low Q T region. This 
also implies that all the above results obtained for the regions (i) and (ii), as well as those 
obtained in §5, are unchanged using the new form (6-13). Note that this is a consequence of 
the fact that (6-13) does not involve any cut-off parameter, like bu m in (6-11), and it is an 
advantage of our approach over the prescription. 

Now we can use (6-13) as the precise form of the resummation formula to investigate its 

*) Actually, we have chosen < 8 < 7r/4 in the above contour C±, anticipating the Gaussian smearing 
factor (9-3) as the corresponding nonperturbative effects. 




(6-14) 
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behavior for Qt < A^cd- Here we mention in particular a remarkable point regarding the 
case Qt ~ 0: In this case, the behavior of (6-13) is controlled by a saddle point in the b 
integration. As shown in §8, the corresponding saddle point is on the real-6 axis with A ~ 1. 
The important role of the saddle point for Qt ~ was pointed out in Refs. 26), 28) and 32), 
using the "old" form, (5-4). Around the saddle point, we have L ~ l/o; s (/i|,) ~ ln(Q 2 / Aq CI) ) . 
Because this implies that all logarithms, L and ln(Q 2 /AQ CD ) ~ l/«s, are regarded as equally 
large for Q ^> Aqcd, the resulting contributions to the resummation formula (6-13) are orga- 
nized in terms of the single small parameter a s , but with a classification of the contributions 
with respect to the order of a s that differs from that of the customary perturbation theory 
that can be used in the region (i), corresponding to Qt > Q- In the Sudakov exponent (5-8), 
one such logarithm, \ii(Q 2 /k 2 ), explicitly appears, and another large logarithm is implicit 
in the integration over k 2 . Suppose that one wants to evaluate the resummation formula 
(6-13) in an approximation of "degree iV," 32 ) meaning that any corrections are suppressed 
by a factor of [1u((5 2 //1qcd)] ~( n+1 \ Then, since two large logarithms multiply the function 
A in (5-8), one must evaluate A in (5-12) to order a^ +2 . Similarly, one needs B to order 
af +1 , C to order af, and, for the running of a s , the (5 function to order af +2 - In particular, 
if one wants an approximate result for (6-13) that will converge to the exact Q T ~ cross 
section as Q — > oo, one needs a degree approximation: A to order a 2 , B to order a], C to 
order a° s , and (3 to order a 2 . 32 - 1 Note that this classification controlling (6-13) for Q T ~ is 
different from that in the customary perturbation theory for Qt > Q, and also from that in 
the all-order resummed perturbation theory, which controls the towers of logarithms (5-1) 
in the kinematical region, Aqcd Qt < Q- This point is discussed in detail in §8. 

Combining the above considerations for the three regions in the resummation formula 
(6-13) with (5-5) and (5-6), we see that this formula represents a complicated object involv- 
ing contributions of various orders and from various distance scales, due to the ubiquitous 
gluons. The impact parameter 6-space formulation allows us to organize these complicated 
contributions into the universal form of (5-5), and this form directly embodies a simple phys- 
ical picture behind the process (5-2): When the transverse momentum Qt of the final-state 
system F is small, the emission of the accompanying radiation is strongly inhibited, so that 
only soft and collinear partons, which have low transverse momenta £>r, are radiated into 
the final state. The quantity da^^iQ 2 ) / d(p represents the hard annihilation of partons, and 
the associated virtual corrections at transverse-momentum scales £>r ~ Q are included in 
the coefficient functions Cji. The form factor e" 5 ^ 6 '^ contains real and virtual contributions 
due to soft [the function Aj(a s ) in (5-8)] and flavor- conserving collinear [the function Bj(a s ) 
in (5-8)] radiation at scales Q > k T > l/b. At very low momentum scales, £>r ^ 1/b, real 
and virtual soft-gluon corrections cancel, because the cross section is IR safe, and only real 



46 



and virtual contributions due to collinear radiation remain, inducing the contributions to 
the coefficient functions Cji(z, a s (bl/b 2 )) and the RG evolution of the parton distributions 

f i/h (x,b 2 /b 2 ). 

§7. The NLL+LO cross section for the tDY 

The (^-differential cross section (5-3) is given by the sum of the resummed compo- 
nent (6-13) and the finite component determined by (5-32). For Q T <C Q, (5-3) is domi- 
nated by the resummed part, da^ TCS ^ / 1 dQ 2 dQ\dyd(j) ) which behaves as ~ l/Q T [ see (5-29), 
(6-5) and (6-6)]. On the other hand, for Q T ~ Q, it approaches the fixed-order result, 
da I dQ 2 dQ^dyd(f)\ fo , since there is no logarithmic enhancement and da^ cs:> /dQ 2 dQ^dyd(f) 
and da^ rcs -y dQ 2 dQ^dyd(f)\ {o almost cancel.*) Therefore, in (5-3) with (5-32), the resummed 
cross section, which is dominant at small Qt, and the fixed-order cross section controlling 
the large Qt region are most naturally matched at intermediate Qt, without double count- 
ing. For example, as discussed in §5.2, the resummed component da^ cs:> jdQ 2 dQ\dyd§ at 
NLL accuracy is exact up to 0{a s {ji 2 R )) when expanded in powers of a s (fi 2 R ). Therefore, the 
NLL resummed cross section should be matched with the fixed-order cross section at LO, 
which is of order a s for Q T > [see the discussion below (5-31)]. We refer to the resulting 
cross section (5-3) as the "NLL+LO" cross section, da NLL+LO j dQ 2 dQ 2 T dyd<\>. Similarly, the 
NNLL resummed cross section is matched with the NLO cross section, giving the result 
for the "NNLL+NLO" cross section, and so forth. Note also that, in order to evaluate the 
NLL+LO cross section with consistent accuracy, the NLO parton distributions have to be 
used as the input for fi/ h (x,fi 2 ). The evaluation of the NNLL+NLO cross section requires 
the NNLO parton distributions. 

7.1. The spin- dependent cross section 

In what follows, we calculate the NLL+LO spin-dependent cross section of the tDY, per- 
forming the matching with the fixed-order result (4-29) obtained in §4. We also "reorganize" 
the resummed component given in (6-13) and (5-5) beyond the CSS form. This is necessary 
for the consistent evaluation of the b integral in (6-13) to the required accuracy. 

Following the general formula described in §§5, 6, we perform the Qt resummation for 
tDY at the NLL level. At this level, we must know the coefficients A {i \ B {1 \ C (1) } 

in (5-12). Then, substituting (5-13), (5-30) and (5-31) into (6-13) and (5-5), the NLL-level 
resummed component of the cross section is given by ^ / cos(20 fel — 0i — 02)^tA nll [see 

*) This cancellation is often very subtle in a numerical sense. See, for example, Ref.54). 



47 



(4-29)], with 

ml dz f 1 dz' 



A T X^(Ql,Q*,y) = [db b -J (bQ T )e s ^ [ ^ I %A T C (z,a s (bl/b 2 )) 

JC 2 Jx° * Jx% * 

xA T C{z\aM/b 2 ))5H^l X f h f^. (7-1) 



Here, the relevant coefficient function is given by 



A T C (z, a s ) =S(l-z) + ^AtC^(z) , (7-2) 



with 



A T C^\z) = C F {^-A ] j 6(1 -z) , 



(7-3) 



and 

S{b,Q) = —^ r -h<®{\)+hW(\) (7-4) 

with (5-23), (5-24) and (5-13) is the exponent (5-22) for the Sudakov quark form factor up 
to NLL accuracy. Above, we have also introduced the shorthand notation for the contour 
integration of (6-13): 

j c db h -UbQ T )W{b- Q,y,cf>)= jf db^H^\bQ T )W(b; Q,y,<j>) 

+ / db h -H^\bQ T )W{b-Q,y,<t>). (7-5) 

J C— 

In (7-1), the parton distributions, as well as the coupling constant a s (b^/b 2 ) associated 
with the coefficient function AtC^\ depend on 6, in addition to the Sudakov factor e s ( b '®\ 
Similarly to the Sudakov exponent S(b, Q) given in (7-4), we can reorganize this b dependence 
in terms of a systematic large logarithmic expansion, using A of (5-21). In order to extract 
this b dependence explicitly, we take the double Mellin moments of (7-1) with respect to the 
DY scaling variables x\ 2 & t fixed Q, 

A T X^%(Q\Ql)= I' dx^xT 1 - 1 f dxl{xl) N ^A T X^{Ql,Q\y) 



JO 



db h -MbQ T )e s ^A T C Nl (a s (b 2 /b 2 )) A T C N . 2 (a s (b 2 /b 2 )) 



x5H Nl , N2 (b 2 /b 2 ) , (7-6) 

where [see (2-4)] 

8H NltN2 (fi 2 ) = / dxix^ 1 ' 1 / dx 2 X2 2 ~ 1 SH(x 1 ,x 2 ; /J 2 ) 
Jo Jo 
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- ^ e 'jfj/h 1 ,N 1 (L J > 2 )fj/h2,N 2 (lJ> 2 ), 

fj/h,N(^ 2 )= [ dxx N - 1 f j/h {x,fi 2 ), (7-7) 
Jo 

with f q/h (x, /j, 2 ) = 5q h (x,/j, 2 ), fg /h (x,/j, 2 ) = 5q h (x,/j, 2 ), and*) 

A T C N (a s ) = / dxx N - l A T C (x, a s ) = 1 + ^ A^. (7-8) 
We rewrite the moment of the coefficient function, A T Cjy (a s (bl/b 2 )), in (7-6) formally as 

A T C N (a s (b 2 /b 2 )) = A T C N (a s (Q 2 )) exp ■ ( 7 ' 9 ) 

Also, note that we obtain the following large logarithmic expansion of the running coupling 
constant (5-19) using (5-21): 45 ) 



In 



a s (Q 2 ) 
<b 2 /b 2 ) 



ln(l - A) + a s (fx R 



Aln(l-A) , _ , (Q 2 \ A 



+ /5 ln 



Po i-A ' V/4/ i-A 



+ 0(« s (4) 2 ) • 

(7-10) 



Substituting (7-8) into the exponent of (7-9), and using (7-10), we find 

h ^g .f)^-A- + oK l^. ( , n) 

Therefore, the exponent on the RHS of (7-9) corresponds to the NNLL level, according to 
the order counting of the contributions in the Sudakov exponent discussed below (5-22). 
Because the NNLL- level term was ignored in the Sudakov exponent (7-4), we can make the 
replacement 

A T C N {a s (b 2 Jb 2 )) - A T C N {a s (Q 2 )) (7-12) 

for the coefficient function in (7-6) at NLL accuracy. 

The b dependence of the parton distributions in (7-6) can be reorganized similarly. In 
this case, the moments of the DGLAP kernel of (5-16) give the corresponding anomalous 



*) In the present case, actually, the moment of the coefficient function, Cji^N, is diagonal in (j, i) and is 
independent of N, as seen from (7-8) and (7-3). However, even when the Cj%,N have non-zero components 
for j 7^ i and depend on 7Y, the results given in (7-11) and (7-12) below hold for each (j, z)-component and 
each moment N of Cji t N- 
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dimensions 7ij,jv(a: s ), and, for the NLO parton distributions relevant to the NLL+LO cross 
section, we have 

^,nM = j\zz N -^^{z,a s ) = ^7giv+ (^)%^, (7-13) 



where the RHS constitutes the perturbative expansion up to the two-loop term, and 7^^ 



and 7^ are given by (A-7) and (A-8) for the transversity distributions. Then, the DGLAP 
equation (5-15) yields the RG evolution, 

fi/hAtf/b 2 ) = J2 U*A%I#, Q 2 )fj/hAQ 2 )i (7-14) 

j=Q,Q 

where Ui^^(b\jb 2 ,Q 2 ) = [1^at(6q/6 2 , Q 2 )]ij denotes the evolution operator, where 

f f b2 ° /b2 dK 2 1 
U N (b 2 /b 2 , Q 2 ) = Pexp I J — lN (a s (K 2 )) \ , (7-15) 

with 7r,-,Ar(a s ) = [j N (a s )]ij, and the symbol U V" denotes the "path ordering" expansion of 
the exponential matrix. Substituting (7-13) into (7-15), the result possesses the well-known 
structure, 57 )' 58) 

U N (b 2 /b 2 , Q 2 ) =[l + (a s (b 2 /b 2 ))] exp | ^ In (^j^y) } [l + O (a s (Q 2 ))] , 

(7-16) 

where the O {a s {b^/b 2 )) term and the O (a s (Q 2 )) term involve the two-loop anomalous di- 
mension matrix "y$ , and their explicit forms for the transversity distributions can be found 
in (A-6). Using (7-14), (7-16) and (7-10) for the parton distributions in (7-6), we find that 
we can replace the parton distributions in (7-6) as 35 - 1 

f l/h ,N(b 2 o/b 2 ) - £ U lJ>N (b 2 /b 2 ,Q 2 )f j/h , N (Q 2 ), (7-17) 

j=Q,Q 

with 

U l3 , N (bl/b 2 ,Q 2 ) = 5 ije R ^ x \ R N (X) = ln(l - A), (7-18) 

up to irrelevant corrections which are down by a s (/xfj) and vanish for A = 0. These corrections 
correspond to the NNLL-level term omitted from the Sudakov exponent (7-4). Now, (7-6) 
reads, to NLL accuracy, 



A T Xxl% 2 (Q 2 , Qt 



1 + -^C F (* 



SH NuN2 (Q 2 )I NuN2 (Q 2 T ,Q 2 ) , (7-19) 



I Ni ,n 2 (Qt,Q 2 ) = ^^Jo(&gT)e s(b ' Q)+ ^ (A)+ ^ (A) . (7-20) 
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In (7-19), the large logarithms L = ln(Q 2 b 2 /6q) and the associated 6-integral is factorized 
into In 1 ,n 2 (QtiQ 2 )> an d therefore, the all-order resummation of the large logarithmic con- 
tributions is now realized at the partonic level. 

If we expand e s ( b >Q)+ R N 1 W+Rn 2 (a) j n p 0wers Q f \ j n the integrand of (7-20), similarly to 
(6-2) and (6-3), the most enhanced logarithmic contributions coincide with the leading tower 
a™L 2n in (6-3), which comes from the expansion of e -( A(1) / 27r ) Q s£ 2 associated with the LL term, 

/ot 8 , of the Sudakov exponent (7-4), and result in the contributions a™ In 2 ™ -1 (Q 2 /Q 2 ^) / Q 2 ^ 
in the cross section, as shown in (6-5). Similarly, noting that [see (7-18)] 

4(i) 5(1) + /' W 0) 
S(b, Q) + R Nl (A) + R N2 (X) = -^-a s L 2 ^ ^« S L + 0(a 2 ), (7-21) 

the next-to-leading tower a™L 2ri_1 , which leads to a™ ln 2n ~ 2 (Q 2 /Qt)/Qt, comes from, e.g., 
the "cross terms" oc (A^a s L 2 ) n ^ 1 x [B^ + 7!?^ + r y^ N ^a s L. For the terms corresponding 
to the next-to-next-to-leading tower, ct™L 2n ~ 2 , which lead to a™ \n 2n ^ 3 (Q 2 /Q^/Q 2 -,, the coef- 
ficients and A^C^ 1 participate, in addition to an< ^ 7^at 2 - As mentioned 
below (6-9), to completely determine the next tower, a™L 2n ~ 4 , the NNLL coefficients A^ 
and -B*- 2 ' associated with /i^(A) of (5-22) and also the coefficient functions and anomalous 
dimensions at the two-loop level, AtC^ n and 7^^, are necessary. 30 ' The NLL resumma- 
tion formula given in (7-19) and (7-20) fully sums up the first three towers of logarithms 
contributing to the cross section, i.e., a™ W n {Q 2 /Q 2 ^) with m — 2n — 1, 2n — 2 and 2n — 3, 
for all n. 

As discussed in §6.1, for Qt > Q, corresponding to the region (i), the resummation is 
irrelevant, and the cross section can be accurately obtained with fixed-order perturbation 
theory. We can switch between the resummed and fixed-order calculations at a certain value 
Qt ~ Q^ but such procedure, introducing ad-hoc boundaries between the large-Q^ and 
small Qi-regions, is not convenient for the actual evaluation of the cross section. For this 
reason, we use (7-19) and (7-20) over the entire region of Qt- The formula (7-20) exponenti- 
ates the functions of A in the integrand, as a result of the large logarithmic expansion, which 
is valid in the regions (ii) and (iii) of §6. Thus, this exponentiated form is effective when 
Qt ^ Q, corresponding to (ii) and (iii), but its behavior cannot be trusted for small b < 1/Q 
{Qt ^ Q)', indeed, L in (5-21) is large for small b as well as for large b, and thus (7-20) also 
exponentiates the unjustified large logarithmic contributions at large Qt- *' This problem 
can be avoided by making the replacement 35 '' 55 ' 

L = In {Q 2 b 2 /bl) -> L = In (Q 2 b 2 /b 2 + l) (7-22) 



*) When we expand the integrand of (7-20) in powers of a s , the enhanced behavior of L for small b is 
harmless, order-by-order in a s [see (5-29), (6-5)]. This is because the b integration in l n of (5-27) converges 
as b — > 0. 
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in the definition (5-21) of A. Through this replacement, the unjustified logarithmic contri- 
butions are suppressed: as b — > 0, we get L — > 0, and e s ( b '^ +RN i ( A ) +Rjv 2 W — > 1. By contrast, 
the resummation of the large logarithmic contributions in (7-19) and (7-20) for b > 1/Q 
is not affected by (7-22), since L and L are equivalent for organizing the large logarithmic 
expansion, as we have L = L + 0(1/ (Q 2 b 2 )). After this replacement, (7-19) and (7-20) read 



a s (Q 2 
2vr 



^H NljN2 (Q 2 )I NljN2 (Q^, Q 2 



L->L 



(7-23) 
(7-24) 



and we denote the double inverse Mellin transform of A T X^^ N JyQ\, Q 2 ) from (iVi, N 2 ) space 
to (xijXjj) space as Z\ T X NLL (Qf,, Q 2 , |/), which represents the new resummed component to 
NLL accuracy. 

The NLL resummed component obtained in this way is now combined with the fixed 
order cross section (4-29). Specifically, the finite component in the NLL+LO cross section is 
determined by (5-32), where we use (4-29) as the first term on the RHS. The result represents 
the NLL+LO cross section of the tDY in the MS scheme, 

ZWa NLL+LO 



dQ 2 dQ^dyd(j) kl 



JV cos (2<j> kl - 0x - <t> 2 ) A T X NLL (Q 2 , Q 2 , y) + A T Y (Q 2 T , Q 2 , y) 



(7-25) 



where 



A T Y(Q 2 T , Q 2 , y) = A T X(Q 2 T , Q 2 , y) + A T Y(Q 2 T , Q 2 , y) - A T X NLL (Q 2 T , Q 2 , y) 



fo 



(7-26) 



Here, (•••)!&> now denotes the expansion up to O (a s (fi 2 R )). Note that A T Y(Q\ 1 Q 2 , y) 
would coincide with AtY(Q^, Q 2 , y), if we did not perform the replacement (7-22) for 
Z\ T X NLL (Q!,Q 2 ,y). We calculate the "subtraction term" Z\ T X NLL | fo in (7-26) explicitly 
[compare with (5-29)] and obtain 



A T X^ LL (Q 2 T ,Q 2 ,y) 



fo 



1 + ^^C F {, 2 - 8) 
2ir 



l Q 5H{xlxl-^ F ) 



+ ^ML Cf ) S H(4,4;ti,)(-Z 2 ^ 3J ;I .) 



2tt 
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+ ; In K 



+ 



1 dz 



o Z 

2 




A T P qq (z) 5H 



where 



db-J (bQ T )\n n (Q 2 b 2 /b 2 + l). 



(7-27) 



(7-28) 



Similarly to the case of X n , (6-14), we can relate this to J c db(b/2)J (bQT){Q 2 b 2 /bl + 
where the contour C can be deformed into that from b = to b = oo along the real 6-axis 
using analytic continuation for rj and Cauchy's theorem. The resulting integral yields [see 
(C-6)] 

J = 5(Q 2 T ), 
b 



11 = -<? 

Zbo 

Q 



x 2 = 



1 u (boQi 
Q- T Kl \-Q 

/v J b Q T ^ \n(Q/Q T ) | 



ir 



(7-29) 
(7-30) 



+00 



Q J Qt b Q 2 T U V Q 
where K n (x) denotes the modified Bessel functions of the second kind, and we have intro- 
duced the notation "[•••] +0O " to indicate the "generalized + distribution" defined between 
< Q\ < 00, as J °° dQ\\- ■ ■ } +00 = (see Appendix B of Ref. 35)). The quantity A^X given 
in (4-31) and (4-32) is identical to (7-27) with X n replaced by I n , and X n is given by X n with 
the replacement L — > L of (7-22). Because the replacement L — > L does not affect the large 
b behavior of the corresponding integrand, the small Qt behavior of I n is the same as that 
of X n . In fact, we can show 35 ) that for Qt Q, we have X n = X n x (1 + 0(Q T /Q 2 )), i.e., 
X n -X n = 0(\n n -\Q 2 T /Q 2 )/Q 2 ) [see (6-4)], and the difference A T X - Z\ T X NLL | fo is given by 
contributions less singular than 1/Q T or S(Q T ) as Q T — > 0.*^ Therefore, AtY differs from 
A T Y only by the 0(a s ) terms, which are less singular than 1/Q T or 5(Q T ) as Q\ — > 0; the 
entire effect of the replacement L — > L is to move a portion of the "finite" component to the 
"singular" component in the cross section. We can rewrite (7-26) as 



A T Y(Q 2 T ,Q 2 ,y) = A T X^\Q 2 T ,Q 2 ,y) 



+ A T Y(Q 2 T ,Q 2 ,y) 



-A T X^(Q 2 T} Q 2 ,y) 



0(a s ), Ql>0 



(7-31) 



*) The fact that the terms proportional to 6(Qj.) cancel out in the difference AtX — Z\T^ NLL |fo is 
actually a nontrivial point. In the present NLL+LO case, they do indeed cancel, as explicitly shown in 
Appendix C. 
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where (• • ■ )\o(a s ) denotes the order a s (fi 2 R ) term in the expansion of (• • • ). Note that A^Y is 
regular over the entire region of Qt, except for a "weak" singularity oc \n(Q 2 / Q T ) / Q 2 , which 
also exists in A^Y .*> The contributions of the first and second terms in (7-31) correspond 
to the LO cross section (4-34), and hence, the NLL+LO cross section (7-25) is actually the 
NLL resummed cross section with the 0(a s ) terms of its expansion subtracted, plus the LO 
cross section. 

One of the advantages of making the replacement L — > L is that the NLL+LO cross 
section defined above satisfies the "unitarity constraint," 35 - 1 which means that (7-25) inte- 
grated over Q T exactly reproduces the Q^-integrated cross section at NLO. Noting that 
e s(b,Q)+R Nl (\)+R N2 (\) = i for 6 = [see (7-24)], we find that the Q T -integral of Z\ T X NLL is 
given by 



/•oo 

/ dQ 2 T A T X NLh (Q 2 T ,Q 2 ,y) 
Jo 



a s (Q 2 j._, 



7T — 



5H{xlxlQ 2 ) (7-32) 



Jo 



2tt 

dQ 2 T A T X(Q T .Q 2 .y). 



'0 

with A T X(Q T: Q 2 : y) given by (4-31) and (4-32). In the second equality, we have used the 
relations, 

i-Q 2 poo 

/ dQ T X n = g n , / dQ T X n = 5 nt o, (7-33) 

Jo Jo 

with g — 1 an d f?i = Q2 — 0, which follow from (6-4) and (7-28), respectively.**) For the 
Z\ T F-term, using (7-26) and (7-33), we find 

dQ 2 T A T Y(Q 2 T ,Q 2 ,y) 

= / dQ 2 T [A T X(Q 2 T , Q 2 , y) + A T Y(Q 2 T , Q 2 , y)] 

Jo 



/•oo 

/ dQ 2 T A T X NLL (Q 2 T ,Q 2 ,y) 
Jo 



fo 



T,max r^T,max 

dQ 2 T A T X(Q 2 T ,Q 2 ,y)+ / dQ 2 T A T Y(Q 2 T , Q 2 , y), 

Jo 



(7-34) 



*) This weak singularity does not matter in da/dQr instead of da jdQ\. 

**) We note that g 3 = — 4£(3), — 0, Q5 — — 48£(5), • • • , and the nonzero contributions in (7-33) for n > 3 
come from the "contact term" proportional to S(Q T ) in (6-4). Therefore, in general cases at NNLL+NLO or 
higher, in principle, we have to take into account those contact-term contributions to satisfy the correspond- 
ing unitarity constraint. To avoid this complication, one convenient method is to perform the matching 
between the resummed and the fixed-order cross sections at Qt > from the beginning. 35 ) 
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where QT,max denotes the maximum value of Qt determined from the partonic kinematics, 

0. - n^ZMMZMZ! (7 . 35) 

^T,max-y • \'™) 

Adding (7-32) and (7-34), we find that the Q T -integral of the NLL+LO cross section (7-25) 
is identical to that of the fixed-order Q T differential cross section (4-29), and hence 



dQ2 A r cfo NLL+LO = A T da 
iT dQ 2 dQ\dyd(j) kl dQ 2 dyd(j) kl ' 



Here, we note that the RHS, i.e. the Qr-integral of (4-29), is counted as the NLO cross sec- 
tion, because the Qr-integrated cross section at LO is of 0(a®), and its partonic subprocess 
is expressed by the tree diagram in Fig. 1. 

7.2. The unpolarized cross section 

The NLL+LO Qr-differential cross section for the unpolarized DY process is obtained 
in the same way as the spin-dependent cross section in §7.1, and the results are given in 
Appendix A of Ref. 56) (see also Refs. 31 )' 32 )). The NLL+LO QT-differential cross section for 
unpolarized DY reads 

d(J NLL + LO 



1JY 



dQ 2 dQ%dyd(j) kl 
where the NLL resummed component X NLL is given by 



X^(Q^Q\y) + Y{Q^Q\y) , (7-37) 



X NLL (Q 2 T ,Q 2 ,y) = I db b -UbQ T )[e s ^ 



2 



H ( x 1 , x 2 ', ^ 2 



<Q 2 ) 

2tt 

+ 



L->L 



(7-38) 



(7-39) 



where the coefficient functions in the MS scheme are given by 

C$}(z) = C F {1 - z ) + C F (^-^j , C$(z) = 2T R z(l - z) , 

with T R = 1/2, and H{x\, x 2 ; ^ 2 ) and Ki :2 (xi, x 2 ; fi 2 ) denote the products of the NLO 
unpolarized parton distributions, 

H(x 1 ,x 2 ;/J 2 ) = ^24 [q hl (x!, /J?)q h2 (x 2 , /J, 2 ) + q hl (x u ^ 2 )q h2 (x 2 , ^m 2 )] , (7-40) 
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K 1 (xi,x 2 ;v 2 ) = J^ e g [QhAxi, fJ?) + q hl (xi,fi 2 )] g h2 (x 2 ,fi 2 ), 



(7-41) 



g 



K 2 (x u x 2 ;^ 2 ) = ^e 2 q g hl (x U [i 2 ) [q h . 2 (x 2 , [i 2 ) + q h2 (x 2 , ^ 2 )] , 



(742) 



where qn(x, fi 2 ) and gh(x, /x 2 ) are the quark and gluon density distributions inside the hadron 
h. In (7-38), the same manipulations as for A T X NLL , such as the replacement (7-22) and 
the contour deformation of (7-5), are applied. Also, the reorganization of the b dependence 
of the parton distributions according to (7-17) and (7-18) should be understood in terms 
of the corresponding anomalous dimensions. Note that the Sudakov exponent S(b, Q) in 
(7-38) for the unpolarized DY process is same as that in (7-1) for the polarized one, given 
by (74). Similarly to A T Y in (7-26), Y in (7-38) is defined using the finite component Y 
of the fixed-order a s cross section. (The explicit expression for Y is given in Appendix A of 
Ref. 56).) Using logic similar to that for (7-32)-(7-36), we can show that (7-37) satisfies the 
unitarity condition 



where the RHS denotes the NLO Qr-integrated cross section for unpolarized DY. 

§8. Asymptotic behavior of the resummed cross section at Q T = 

In §6, we discussed how the resummation formula given in (54) and (6-13) is controlled 
in each of three regions (i), (ii), and (iii) of the impact parameter b, corresponding to the 
kinematical regions, Qt~Q, ^Iqcd *C Qt~ Q and Qt~^iqcd, respectively. In the region 
(iii), the Sudakov form factor e s( - b '^ with (5-22) cannot be expanded in terms of A, since 
|A| = 0{1). In this section we show that the behavior of (6-13) in the region (iii), in particular 
its Qt — > limit, is controlled by the saddle point in the b integration. Such analysis has 
been done in Refs. 26) and 32) for the LL resummation formula. Here we extend the analysis 
to the case of the NLL-level resummation formula with (7-23) and (7-24). It is possible to 
carry out the corresponding extension on the basis of the present formalism, which realizes 
the resummation at the partonic level, as discussed in §7. 

Let us start from the NLL resummation formula (7-24) with the exponentiation of the 
functions of A in the integrand. For simplicity, we fix the renormalization scale as /ir = Q 
in the following. When Q is large enough that a s (Q 2 ) <C 1, the b integral in (7-24) at 
small Qt is dominated by a saddle point, which is determined mainly by the LL term in 
the Sudakov exponent (74), h^(X)/a s (Q 2 ) — > oo. In this case, the contributions to the b 
integration from very short (|6| <C l/Q) as well as very long distances (|6| ^> 1/tIqcd) along 
the integration contour C± are exponentially suppressed [see (643), (7-5)]. This allows us 




(743) 
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to give up the replacement (7-22) in (7-24). Also, we may ignore the integration along the 
two branches, b = b c + e ±ie t with t G {0, oo} in C±, when b c is sufficiently large, but less 
than the position of the singularity in the Sudakov exponent, bi = (bo/Q)e 1 ^ 2f3 °°' a ^ 2 ^ . (In 
fact, the relevant integrand has a nice saddle point well below b^ (above 0) for Q ^> Aqqd, 
as seen from (8-6) and (8-11) below.) Then (7-24) reads, up to the exponentially suppressed 
corrections, 



f bc b 

In u n 2 (Qt,Q 2 ) = / db -J (bQ T ) exp 
Jo * 



lf^ + h^(X)+R Nl (X) + R N2 (X) 



i) 



with b c < b L . To evaluate the integral in (8T) with the saddle-point method, we consider 
the simplest case, with Q T = 0, and change the integration variable from b to A of (5-21), as 

WQt = 0, Q 2 ) = An2 f° (n2 , f C dAe-C (0) (W (1) W+*-i W+^W, (8-2) 

where A c = (3 a s {Q 2 ) ln{Q 2 b 2 c /b 2 ) and 

AQ) (X)= * fe(°)(A) 

C (A) - Poa s (Q 2 ) a s (Q 2 Y ^ 

Here, C^(^) denotes the leading term oc l/a s (Q 2 ) in the exponent of (8-2), while the re- 
maining terms, ft/ 1 ) (A) + Rn 1 (X) + Rn 2 (X), collect the contributions of the NLL level, which 
are down by a s (Q 2 ) from the leading term in the relevant region, A ~ 1. The precise position 
of the saddle point in the integral of (8-2) is determined by the condition 

-C (0) '(A) + h^'(X) + R Nl '(X) + R N2 '(X) = 0. (8-4) 

Now we express its solution as A = Xsp + AXsp, where Xsp is the solution of (^'(X) = 0; 
i.e., Asp satisfies 

1 " S-T^T- = < 8 ' 5 > 

ZTTfjQ 1 - X S p 

and AXsp denotes the shift of the saddle point due to the subleading terms in (8-4), h^'(X) + 
R Nl '(X) + R N2 '{X). From (8-5), we get 

Xsp = ttt (8-6) 

2ir[3 + 4 } 

which is independent of Q and coincides with the saddle point at the LL level discussed 
in Refs. 26), 32) and 51). Substituting A [ g ] = 2C F from (5-13), we obtain X SP ~ 0.6, 
which is 0(1) as it should be. Contrastingly, the shift, AX S p = [h^'(X S p) + Rni(^sp) + 
Rn 2 '(^sp)]/(^"(^sp) + •••) is 0(a s (Q 2 )) or higher. Up to the exponentially suppressed 
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corrections, the saddle point evaluation of (8-2) around A = Xsp + AXsp = uo can be 
performed in the usual way. Defining £(A) = (^(X) — h^(X) — R^X) — Rn 2 (X), we obtain 



^Vi,AT 2 (0, Q 2 ) 



1.2 rco 

" / dAe- rH/2](A - w)2 



4Q 2 /5 a s (Q 2 ) 

x (l - 2^"»(A - c) 4 + l[f ' (c)] 2 (A - c) 6 + ■ 



4Q 2 /? a s (Q 2 ) V Z» V § [£"M] 2 24 [{"(a;)] 3 



(8-7) 



Because ~ l/a s (Q 2 ) for a s (Q 2 ) <C 1, the second and third terms in the parentheses 
on the RHS of (8-7) are of 0(a s (Q 2 )) ) and similarly it is straightforward to show that the 
subsequent terms denoted by the ellipses are of O (a s (Q 2 ) 2 ) or higher. We ignore the second 
and following terms in the parentheses, because they correspond to the NNLL or higher-level 
contributions that were omitted from our starting point (7-24). As discussed below (5-22), the 
NNLL contributions are associated with the terms a" +1 L n oc a s X n (n > 1) in the exponent 
of the resummation formula. More precisely, in the regions (i) and (ii) of §6, those NNLL 
contributions behave as ~ a 2 L, a^L 2 , . . ., while at the saddle point A = uo = X S p + AX S p ~ 1, 
they would behave as ~ a s . A similar situation occurs for the NNNLL and higher-level 
contributions. As mentioned in §6.3, this indicates that, for Qt ~ in the region (iii), all 
contributions are completely controlled by the single small parameter a s (Q 2 ), and we have 
to use a classification with respect to the order counting of the terms that is different from 
the customary and all-order resummed perturbation theory used in the regions (i) and (ii). 
In particular, when we start with the NLL-level resummation formula and go to the Qt ~ 
region, we should ignore the contributions that are of the same order as the 0(a s ) terms in 
the exponent of the corresponding formula. 32 )> 56 ) 

Substituting uo = Xsp + AXsp and £ = — — R^ — R^ 2 into (8-7), and expanding 
the result in terms of AXsp, we find 



f AT ... (ft Q2) - b 271 r -\C m ^sp)+h( 1 HX 3P )+R N A^sp)+RN,(X S p)} (O.Q) 

to NLL accuracy. Note that using C^'(Xsp) = 0, the shift AX S p yields only NNLL-level 
corrections. We substitute this into (7-23) and combine the result with the coefficient func- 
tion and the parton distributions. Then, using the same logic as above, the 0(a s (Q 2 )) term 
in the coefficient function in (7-23) is also ignored as [1 + a s (Q 2 )CF{^ 2 — 8)/2n] — > 1, and 
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thus we have 

A T X^ L N2 (0,Q 2 ) = 5H NuN2 (Q 2 )I Nu nMQ 2 ) • (8-9) 
Then, performing the double inverse Mellin transformation to the (x±, x^) space, we obtain 56 ) 

A T X NLL (0,Q 2 ,y) 

{4Q*(3oa s m\j C(°)"(A S p) e ) V x ' 2 ' b 2 SP J ' 

(8-10) 

with 

bsp = ^s P /[2^s(Q 2 )] j ( 8 .H) 

Q 

which is the value of b corresponding to \sp in (8-6). Here, we have used the fact that 
qRn^sp) anc [ 6 Rn 2 (*sp) j n (g.g) can j-jg identified with the NLO evolution operators from the 
scale Q to b /bsp, to the present accuracy. This yields 5H(x\, x%; Q 2 ) with the scale shifted 
8,sQ^b /bsP [see (7-17), (7-18)]. 

The saddle-point evaluation of the NLL resummation formula (7-38) for the unpolarized 
case can be performed similarly, and the result is given by 56 ) 

l NLL (0,Q 2 ,y) 

- ( h l / 2lT r -^H\sp)+hW(X SP )\ rr ( r 0. 3_\ 

\±Q 2 h* s {Q 2 )\i C^iXsp) J V 2 ' b 2 SP ) ' 

(8-12) 

up to NLL accuracy, which corresponds to the above result (8-10) with the replacement 
5H(xi,X2', bl/b 2 SP ) — > H(xi, x®; bl/b 2 SP ). At the present accuracy, the gluon distribution 
decouples for Q T 0, since the corresponding coefficient function Cqg in (7-39) is of order 
a s (Q 2 )- 

We emphasize that the results (8T0)-(8-12) are exact, up to the NNLL corrections corre- 
sponding to the 0(a s (Q 2 )) effects. We also see that these results give the explicit realization 
of the "degree 0" approximation discussed in §6.3 for polarized and unpolarized DY pro- 
cesses. The common factor for both the polarized and unpolarized results, given by the 
contribution in the parentheses of (8-10) and (8-12), represents "very large perturbative ef- 
fects" due to the universal Sudakov factor. Substituting (8-6) and the running coupling 
constant at two- loop level [see (5-19)] into this common factor, we obtain, to the present 
accuracy, 

6 / 271 r -C^(Xsp)+h^(\sp) 

4Q 2 /WQ 2 )V C<°>"(Asp) 
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2v^(a + l) ^qcd V / 

x [ln(Q 2 /4cD)] (/3lM ' )[1 " aln(1+1/a)1 e fe(1)(A ^ ) , (8-13) 

with a = A^/(2ty(3 ) ~ 0.6. Here, the factors in the second line reproduce the well-known 
asymptotic behavior for Q ^> Aqqb, derived from the conventional LL- level saddle-point 
calculation 26 ) (see also Refs. 51) and 52)), while those in the third line represent the NLL- 
level effects. We note that the subleading, NLL contributions do not disappear in the 
asymptotic limit Q — > oo; e.g., e /l(1>( - Asp - ) provides the 0(1) constant factor [see (5-24)]. In 
fact, according to the discussion of the "degree N" approximation in §6.3, our results (8-10)- 
(8-13), representing the degree approximation, are exact in the asymptotic limit, while the 
"degree —1" approximation, which corresponds to the LL resummation, is not. 32 ) 

In our saddle-point formulae (8-10) and (8-12), the "new scale" b^/bsp participating in the 
parton distributions is another NLL-level effect, which does not disappear in the asymptotic 
limit Q — > oo. In fact, for the polarized case (8-10), the corresponding contribution is 
associated with the 0(1) terms, (Xsp) + Rn 2 {^sp), in the exponent of (8-8). Substituting 
(8-6), we can express (8-11) as 

7^ =^QCD (-r-)^ 1 [ln(g 2 M 2 Q CD)]"^, (8-14) 
Osp \ a qcd/ 

where the first two factors on the RHS coincide with the corresponding LL-level expres- 
sion 32 )' 51 ) Q f ho/bsp, while the third one, involving the logarithms of Q 2 / Aq CD , represents 
the NLL-level modification induced by the two-loop running of the coupling constant in 
(8-11). Because a ~ 0.6, this implies bo/bsp ~ Aqqd (Q / ^QCD) a ^ (a+1 ' ) — a/^qcdQ) modulo 
logarithms. Therefore, our saddle-point formulae (8-10) and (8-12) are directly applicable to 
the production of very high mass DY pairs such that bo/bsp ^> Aqcd- In fact, in such cases, 
the integral over b in our resummation formula (7-24) is actually dominated by the region 
of b near b$p (<^C 1//1qcd)), where the form of the integrand in (7-24) is accurate [see the 
discussion below (5-10)]. Practically, however, the polarized DY experiments at RHIC, J- 
PARC, GSI, etc., will measure dileptons with masses Q ~ several GeV, where bo/bsp may be 
of order 1 GeV, and our formulae (8-10) and (8-12) are not useful for analyzing those cases. 
Apparently, a similar problem arises at the LL-level, 26 )' 32 )' 52 ) and the present result indicates 
that the higher-order perturbative corrections at NLL level cannot solve this problem. 

The root of this frustrating conclusion is in the fact that the integrand of our resummation 
formula (7-24) is composed of purely perturbative quantities and is inaccurate in the long- 
distance region (iii), |6| > Aq CD . We have to introduce nonperturbative effects relevant to 
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such a large b region in order to treat the small Qt behavior in the tDY for the kinematics 
of the relevant experiments. We discuss this point and the application of the results to 
phenomenology in the next section. 

§9. Predictions for the transversely polarized Drell-Yan process at RHIC and 
J-PARC in the NLL+LO Q T -resummation framework 

On the basis of the QCD factorization and the transverse- momentum resummation frame- 
work, we have obtained all the necessary analytic formulae of the higher-order perturbative 
corrections for the QCD prediction of the dilepton Qt spectrum and the spin asymmetry 
in the transversely polarized DY process. In this section, we specify the nonperturbative 
inputs necessary for estimating them and calculate those observables numerically with the 
kinematics of an ongoing experiment at RHIC and a possible future experiment at J-PARC. 
A detailed numerical study has already been given in Ref. 56), and here we highlight the im- 
portant results to show how the soft-gluon resummation contributions control the behavior 
of the Qt spectrum and the double transverse-spin asymmetry, especially in the small Qt 
region. We also analyze the double-spin asymmetries, utilizing the saddle-point evaluation 
of the resummation formula. 

9.1. The Qt- differential cross sections 

First, we calculate the Qr-differential cross sections (7-25) and (7-37) for polarized and 
unpolarized Drell-Yan processes, respectively, with the kinematics of pp collisions at RHIC. 
For simplicity, here and below we denote the parton distributions of the proton as [see (2-4), 

(A-l), (740)-(742)] 

5q(x, /i 2 ) = 5q h=p (x, /i 2 ), q(x, /i 2 ) = q h=P (x, fi 2 ), g(x, /i 2 ) = g h=p (x, /z 2 ), (9-1) 

omitting the subscript, and also we denote the combination of the azimuthal angles appearing 
in (7-25) as 2<p kl - fa - fa = 20. 

To calculate the unpolarized cross section (7-37), we need nonperturbative inputs for the 
unpolarized quark and gluon distributions, q(x,fi 2 ) and g(x,fi 2 ); we use the NLO GRV98 
distributions for them. 64 ) To calculate the polarized cross section (7-25), we need nonper- 
turbative inputs for the transversity 5q(x,/j, 2 ). Recently, Anselmino et al. 60 ^ gave the first 
estimation of the transversity distributions from semi-inclusive DIS (SIDIS) data, combined 
with e + e~ data for the associated (Collins) fragmentation function. They performed a global 
fit of the SIDIS spin-asymmetry data with the corresponding LO QCD formula, assuming 
that the antiquark distributions vanish, i.e. d~q(x,fi 2 ) = 0, and obtained the LO quark 
transversity distribution. Thus, their results cannot be used for the present calculation of 
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tDY at the NLL+LO accuracy, in which the NLO transversity distributions for both the 
quark and antiquark are necessary. At the moment, no other experimental information on 
the transversity exists*- 1 Theoretically, a nontrivial relation among the transversity 5q(x), 
unpolarized q(x), and longitudinally polarized Aq(x) distributions is provided by the Soffer 
inequality, 63 ^ 



for each flavor of quark and antiquark. In this paper, we use a model of the NLO transversity 
distributions which satisfy the Soffer inequality, following Refs.20) and 22). In this model, 
we saturate the Soffer bound of (9-2) at the low input scale fio — 0.6 GeV, as 5q(x,fj,o) = 
[q(x, Ho) + Aq(x, /Jq)]/2. Here, we cannot determine the signs of the distributions from 
(9-2), and in this calculation, we have chosen all the signs to be positive. For the input 
functions on the RHS of this formula, q(x, /Jq) and Aq(x, /1q), we use the NLO GRV98 64 ) and 
GRSV2000 ("standard scenario") 65 -* distributions, respectively. Then we evolve 5q(x,Ho) to 
higher scales, according to the NLO DGLAP kernel for the transversity as (A-6). 

We also introduce another nonperturbative contribution to the QT-differential cross sec- 
tions. As mentioned in §8, in the kinematical regions for most experiments, including those 
at RHIC and J-PARC, the small Qt behavior of our resummation formula (7-24) is domi- 
nated by the contribution from \b\ > 1/Mqcd i n the ^ integral, which corresponds to the region 
(iii) of §6.3. However, the integrand of (7-24) involving purely perturbative quantities is not 
accurate for such a large |6| region, and the corresponding long-distance behavior has to be 
complemented by relevant nonperturbative effects. Formally, those nonperturbative effects 
play a role of compensating for the ambiguity that the prescription for the b integration 
in (7-24) needed to avoid the singularity in the Sudakov exponent S(b, Q) of (7-4) is not 
unique, such as the b* prescription of (6-10), (6-11), and the present prescription with (6-13). 
Therefore, following Refs. 32), 45) and 35), in (7-24) we make the replacement 



with a nonperturbative parameter g NP . Because exactly the same Sudakov factor, e s ( b '®\ 
participates in the corresponding formula for the unpolarized case, (7-38), as noted in §7.2, 
we perform the replacement (9-3) with the same nonperturbative parameter g^p in the 
NLL+LO unpolarized differential cross section. This may be interpreted as assuming the 
same "intrinsic transverse momentum" of partons inside nucleons for both polarized and 

*) For a recent result of a lattice QCD simulation of the first and second moments of the it- and d-quark 
transversity distributions, see Ref. 61). The transversity distributions fitted to the prediction of the chiral 
quark soliton model of nucleons are given in Ref. 62). 




(9-2) 



,S(b,Q) p S(b,Q)-g NP b 2 



(9-3) 
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unpolarized cases, corresponding to the Gaussian smearing factor of (9-3). For the present 
calculations, we use g^p = 0.5 GeV 2 , which is consistent with a study of the Qt spectrum 
in the unpolarized case. 36 - 1 

Assuming e~ 9Npb2 — > 1 in customary perturbation theory, it is straightforward to show 
that all results derived in §7 are unchanged under the replacement (9-3). In particular, 
the NLL+LO QT-differential cross sections (7-25) and (7-37) with (9-3) satisfy the unitarity 
condition, given in (7-36) and (7-43), exactly. 

For all the following numerical calculations, we choose = for the azimuthal angle of 
one of the leptons, //p = /ir = Q for the factorization and renormalization scales, and b c = 
and 9 = ^7r for the integration contour C± defined below (6-13). 




Fig. 2. The spin-dependent and spin- averaged differential cross sections for tDY: (a) 
A^da /dQ 2 dQT<iyd4> and (b) da /dQ 2 dQTdyd(p, as functions of Qt with RHIC kinematics, 
y/S = 200 GeV, Q = 5 GeV, y = 2 and <j) = 0, with g NP = 0.5 GeV 2 . 

The tDY (^-differential cross sections for the spin-dependent and spin-averaged cases 
are shown in Fig. 2(a) and Fig. 2(b), respectively, for \/S = 200 GeV, Q = 5 GeV and 
y = 2, corresponding to the detection of dileptons with the PHENIX detector at RHIC. 
The solid curve in Fig. 2(a) shows the NLL+LO differential cross section (7-25), multiplied 
by 2Qt- The dot-dashed curve represents the contribution from the NLL resummed com- 
ponent Z\ T X NLL in (7-25).*) The LO cross section (4-34) is depicted by the dashed curve, 
which becomes large and diverges as Qt — > 0, due to the singular terms a s \n(Q 2 /Q T )/Q T 

*' In principle, we should take the "degree 0" approximation for the numerical calculations when Qt ~ 0, 
as discussed in §6 and §8. But we do not make the corresponding replacement for the coefficient functions, 
Cji{z, a s ) — > SjiS(l — z), for Qt ~ in our numerical calculations, because the effect of the replacement is 
invisible in the solid as well as dot-dashed curve for Qt ~ in Fig. 2(a), and also in Fig. 2(b) below. For 
the spin asymmetry, which is defined as the ratio of the cross sections in Figs. 2(a) and (b) and is discussed 
in §9.2, the replacement would increase the NLL+LO asymmetry (9-4) and the NLL asymmetry (9-5) for 
Qt ~ in Figs. 3-5 by about 5%. 



63 



and a s /Q T in AtX^ of (4-32). This singular behavior disappears in the NLL resummed 
component AtX nll when we sum up the recoil logarithms to all orders in a s . The NLL 
resummed component dominates the NLL+LO cross section in the peak region of the solid 
curve, i.e., at intermediate as well as small values of Qt- As a result, the NLL+LO cross 
section (7-25) is finite and well behaved over all regions of Qt- We also show the LL result, 
represented by the two-dot-dashed curve, obtained from the corresponding NLL result by 
omitting the contributions from the NLL terms, i.e., h^-'(X), Rn 1 (X), Rn 2 (X) in (7-24) and 
o. s {Q 2 )Cf{^ 2 — 8)/27r in (7-23). The LL contributions are sufficient to obtain a finite cross 
section, suppressing the cross sections for Q T « 0. 

The curves in Fig. 2(b) are the same as those in Fig. 2(a), but for the unpolarized cross 
sections; e.g., the solid curve is obtained by calculating (7-37), and the dot-dashed curve is 
the contribution from the NLL resummed component X NLL in (7-37). The behavior of each 
curve is similar to that in the polarized case. 

Comparing the NLL results with the corresponding LL results in Figs. 2(a) and (b), we 
find that the contributions from the NLL terms considerably enhance the cross sections in 
both polarized and unpolarized cases, while the enhancement is larger for the former case: 
The "universal" term h^'(X) gives a similar (enhancement) effect for both channels, while 
the other NLL contributions associated with the evolution operators and the 0(a s ) terms 
in the coefficient functions give different effects in the polarized and unpolarized cases [see 
(7-18), (7-3) and (7-39)]. 




1 1.5 2 
Q T (GeV) 

Fig. 3. The asymmetries for RHIC kinematics, y/S = 200 GeV, Q 
obtained from each curve in Fig. 2. 



5 GeV, y = 2 and 
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9.2. The double transverse- spin asymmetries 

We consider the double transverse-spin asymmetry, defined by the ratio of the NLL+LO 
cross sections (7-25) and (7-37), 

w TT (Qt) - ■= cos(20J v-ntt/^2 ^2 \ , <>ira r,i \ ' ^ > 

2 X NLL (Q T ,Q 2 ,y) + Y(Q T ,Q 2 ,y) 

as a function of Q T - Figure 3 plots the asymmetries for the tDY at RHIC, given by the ratio 
of the value of each line in Fig. 2(a) to that of the corresponding line in Fig. 2(b); that is, 
the solid curve gives the NLL+LO asymmetry £^tt L+L °(Qt) of (9-4), and the dot-dashed 
curve gives the NLL result, defined as 

< T LL (Q T ) = icos(20) ^ XNLL(g ^ g2 ^ ) . (9-5) 

The NLL+LO asymmetry £^^ h+h ° (Qt) in Fig. 3 is almost flat in the peak region of the 
NLL+LO cross sections in Fig. 2. This flat behavior is dominated by the NLL resummed 
components, and it reflects the fact that the soft gluon effects resummed into the Sudakov 
factor e 5 ^'°) with (7-4) are universal to the NLL accuracy between the numerator and the 
denominator of s^ L (Q T ) of (9-5). * } A slight increase of ^ T LL+LO (Q T ) as Q T -> is due 
to the terms oc \n(Q 2 /Q T ) in the "regular components" A T Y and Y in (9-4). The dashed 
curve in Fig. 3 shows the LO asymmetry ^tt(Qt), defined by the ratio of (4-34) to the 
corresponding LO cross section for unpolarized DY process, 

, LO ,„ . 1 ,„ ^tA-W(Q?, Q\ y)\ Ql>a + A T Y[Ql, Q\ y) 

(Qt) = 2«»™ PlWk^iW • (9 - 6) 

where X^ denotes the singular part of the 0(a s ) unpolarized cross section, 31 ) defined simi- 
larly to A T X^ in (4-32). Although both the numerator and denominator of (9-6) diverge as 
Qt — ► 0, as shown in Figs. 2(a) and (b), £^tt(Qt), their ratio, appears to be finite. However, 
this LO result decreases as Qt increases, and it is much smaller than the NLL+LO result, 
indicating that the soft gluon resummation is crucial for the prediction of the asymmetries 
at small Qt and intermediate Qt- The two-dot-dashed curve in Fig. 3 shows the LL result 
[see (7.23), (7.24), (7.38)], 



*) We note that the shape of the Qt spectra in Figs. 2(a) and (b) is actually sensitive to the value 
of the nonperturbative parameter gxp of (9-3) (see Ref. 23)). But the corresponding g^p dependence 
almost cancels between the numerator and denominator in the asymmetries of (9-4) and (9-5) in the range 
g NP = 0.3-0.8 GeV 2 . 56 ) 
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which is constant in Qt, because the soft-gluon corrections in the numerator and denom- 
inator exactly cancel at the LL level [see also (6-6)]. As expected from Fig. 2, £/tt L (Qt) 
is larger than £&tt(Qt), due to the effects of the NLL contributions. In particular, among 
the NLL-level effects noted at the end of §9.1, the contributions associated with the evolu- 
tion operators play an important role to amplify £^tt L (Qt) i n comparison with £/tt(Qt)'- 
In the resummation formula (7-24), the large logarithmic contributions due to the Su- 
dakov factor in (7-4) and (5-8) strongly enhance the contributions from the region b ~ 
1/Qt in the b integration, where the evolution operator (7-18) causes the scale of the 
parton distributions to be b /b ~ Qt, instead of Q. Therefore, the numerator of (9-5) 
is dominated by the transversity distributions at the scale ~ Qt- A similar mechanism 
arises in the denominator, given by (7-38), so that, after the cancellation of the univer- 
sal Sudakov factor between the numerator and denominator, the dominant contribution 
to (9-5) is given by ^ T LL (Q T ) ~ [cos(20)/2] 5H(x\, x° 2 ; Q T )/H(xl, x° 2 ; Q 2 T ). Note that the 
transversity distributions in (2-4) obey QCD evolution of the non-singlet type, described 
by (A-6). Thus, as the scale jj, increases, the antiquark transversity distribution 5q(x,fi 2 ) 
in the small x region grows less rapidly than the unpolarized sea distribution q{x,ji 2 ) in 
(7-40). The latter grows rapidly for small x through mixing with the gluon distribu- 
tion. Therefore, the ratio of the sea distributions 5q(x, fi 2 )/q(x, ji 2 ) at small x becomes 
enhanced as the scale \i is reduced from Q to Qt- This explains why the NLL asymmetry 
^tt L (Qt) ~ [cos(20)/2] 8H(x\, x 2 ; Q T )/H(x° l) x 2 ; Q T ) is enhanced in comparison with the 
LL result £^tt(Qt) appearing in (9-7) in the present case, with = 0.185 and x 2 = 0.003 
in (4-7) [see also the discussion in §9.3 using the saddle-point formula (9T3) below]. 

Now we compare these asymmetries in Fig. 3 with the asymmetry for the Qr-mtegrated 
cross sections A T t, defined as 

JdQ 2 T [A T da NLL+LO /dQ 2 dQ 2 T dyd(j)] 
Att= j dQ 2 T [da NLL + LO /dQ 2 dQ 2 T dyd<j)} 

_ [A T da/dQ 2 dyd<P] _ 1 dHjx ^ x° 2 ; Q 2 ) + ■■■ 

[da/dQ 2 dyd<P\ 2 C ° Sl 9) H(a%, x° 2 ; Q 2 ) + • • • ' 1 ' 

where we have used the fact that our NLL+LO QT-differential cross sections satisfy the 
unitarity conditions (7-36) and (7-43) exactly, and the ellipses on the RHS stand for the NLO 
[(9(a s )] correction terms in the MS scheme. Therefore, Att defined in this way coincides 
with the NLO asymmetry calculated in Ref. 20).*) Using the present nonperturbative inputs, 
we get A T t = 4.0%. This is close to the value of ^tt(Qt) an d is smaller than that of 

*) In Ref. 20), the corresponding asymmetry is denned through a certain integration over <fi, and it is 
equal to (9-8) with the formal replacement cos(2</>) — > 2/tt. Also, the input parton distributions in Ref. 20) 
are not completely the same as ours (see the discussion in Ref. 56)). 
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^4 N T LL+LO (Q T ) in the flat region, which is reasonable, because Att in (9-8) coincides with 
£?tt(Qt) in (9-7), up to NLO QCD corrections. 
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Fig. 4. The NLL+LO J3^ L+LO (Q T ) m (94) with (9-3) using g NP = 0.5GeV 2 with RHIC kine- 
matics: (a) VS = 200 GeV, y = 2 and <f> = 0; (b) VS = 500 GeV, y = 2 and <j> = 0. 



Figures 4(a) and (b) show the NLL+LO asymmetries s^^ +lj0 (Q T ) with RHIC kine- 
matics, \/S = 200 GeV and v^S" = 500 GeV, respectively, with various values of Q and y = 2. 
The dashed curve in Fig. 4(a) is the same as the solid curve in Fig. 3, and similarly to this 
dashed curve, all other curves in Figs. 4(a) and (b) exhibit the characteristic flat behavior 
in the small Qt region. We mention that the dependence of £/^ L+LO (Q T ) on Q comes 
from the small- x behavior of the relevant parton distributions: A smaller Q corresponds 
to a smaller x\ 2 — e ±y Q/\/S, so that the small-x rise of the unpolarized sea distributions 
enhances the denominator of (9-4). A similar mechanism also explains why the results for 
v^S* = 500 GeV are smaller than those for y/S = 200 GeV. Comparing with the NLO Att 
in (9-8), we fmd 56) that the values of £^ l+lo (Qt) in the flat region of Figs. 4(a) and (b) 
are larger by 20-30% than the corresponding values of Att- 

Next, we present the results for J-PARC kinematics, i.e., the results for the tDY in 
pp collisions at the moderate CM energy y/S = 10 GeV. Fi gure 5(a) shows the double 
transverse-spin asymmetries for v^S* = 10 GeV, Q = 2 GeV and y = 0, and the curves have 
the same meanings as the corresponding curves in Fig. 3. The results for J-PARC have 
features similar to those in Fig. 3 for RHIC: The NLL+LO asymmetry £/^ L+LO (Q T ) is flat 
in the small Qt region, and it is dominated by the NLL asymmetry £^tt L (Qt)] the NLL 
asymmetry is enhanced in comparison with the LL result, and the fixed-order LO result is far 
below the others. We note that £^ L+LO (Q T ) ~ 15% in the flat region is again larger than 
the corresponding prediction using (9-8), Att — 12.8%. Figure 5(b) shows the NLL+LO 
asymmetries ^^ L+LO (Q T ) with J-PARC kinematics, v 7 ^ = 10 GeV, Q = 2, 2.5, 3.5 GeV, 
and y = 0, and the solid curve is the same as the solid curve in Fig. 5(a). We observe that 



67 



30 

g 

^20 
© 

< 10 



1 

(a) 


| , | 

NLL+LO 




NLL 




■ - ■ LL 




-- LO 




X 




N. 

"S. 

1 





0.5 1 1.5 

Q T (GeV) 



30 



&20 



10- 



1 

_ (b) 


1 , 1 


Q = 2GeV 




-- Q = 2.5GeV 




■-■ Q = 3.5GeV 






1 





0.5 1 1.5 

Q T (GeV) 



Fig. 5. The asymmetries with J-PARC kinematics, (a) The asymmetries obtained similarly to 
Fig. 3(a), with = 10 GeV, Q = 2 GeV, y = and <j> = 0. (b) The NLL+LO ^ L+h °{Q T ) 
with VS = 10 GeV, y = and <p = 0. 



the asymmetries stf TT + (Qt) are flat and are approximately 15% for all Q; we note that 
this value of £/^ L+LO (Q T ) is larger by 20-30% than the corresponding values of Att in 
(9-8). 56 ) Moreover, the asymmetries found at J-PARC are in general larger than those at 
RHIC energies, because the parton distributions for larger x® 2 >0.2 are probed, where the 
enhancement mechanism of the denominator of (9-4) by the unpolarized sea distributions, 
as noted in Figs. 4(a) and (b), is absent.*) 

Finally, we emphasize that in all cases considered in Figs. 3-5, we have observed that 
,< T LL+LO (Q T ) in (9-4) in the flat region is considerably larger than the corresponding Att 
in (9-8). Indeed, this is generally true for the tDY in pp collisions, regardless of the specific 
kinematics or the detailed behavior of nonperturbative inputs, because this phenomenon is 
mainly governed by the partonic mechanism associated with the soft gluon resummation at 
the NLL level, as demonstrated in Fig. 2. By contrast, apparently the absolute magnitudes 
of both j4t L+L °(<2t) and A TT are influenced by the detailed behavior of the input parton 
distributions, in particular, at RHIC energies, by their small-x behavior. 

9.3. The saddle-point formula incorporating nonperturbative smearing 

In the numerical results presented in §9.2, we have observed the universal flat behavior 
of £^tt L+L °(Qt) at small Qt, where the bulk of the DY pairs is produced. We have also ob- 
served that in this flat region, the numerator and denominator of (9-4) are dominated by the 
NLL resummed components, AtX nll and X NLL , respectively, and as a result, i^^ L+L0 \Qt) 
is well approximated by &^tt L (Qt)- This suggests that the experimental value of the double 

*) It is expected that even larger asymmetries will be observed in future pp experiments at GSI. 67 )' 68 ) 
The application of our Qt resummation formalism to pp collisions will be presented elsewhere. 69 ' 



68 



transverse-spin asymmetry obtained from the data in the flat region can be compared with 
the theoretical value of (9-5) directly. Noting that the flat behavior of the NLL asymmetry 
(9-5) can be extrapolated to Qt — > 0, as demonstrated in Figs. 3 and 5(a), we can use the 
relations 

^tt(Qt) - < t LL (Qt) - < T LL (0) = - cos(20) T ~ (9-9) 

for Q T in the flat region, where £^tt(Qt) denotes the experimental value of the double 
transverse-spin asymmetry. In §8, the Qt — * limit of the NLL resummed components, 
A T X NLL (0, Q 2 , y) and X NLL (0, Q 2 : y), has been evaluated analytically using the saddle-point 
method for the case without the replacement (9-3). In view of (9-9), it would be interesting 
to extend the results in §8, incorporating the nonperturbative smearing factor as in (9-3). 56 ) 
Starting from (7-24) with the replacement (9-3), we can use logic similar to that used in 
§8 in order to obtain (8-1), which now reads, at Q T = 0, 



l N u N 2 



(0,Q 2 ) = \ ! dbbexp 
* Jo 



MM + ftW(A) + R N (A) + R N2 (X) - g NP b 2 



II" ■ > A ) 

(9-10) 



Because b 2 = (b 2 ) /Q 2 )e x l <y P° aa<yQ2 ^ using (5-21), the nonperturbative function produces a 
"double exponential" behavior, and the saddle-point evaluation with such an integrand is a 
nontrivial task. Nevertheless, we can check numerically that the integrand of (9-10) has a 
nice saddle point well below b c (above 0) for the kinematics of interest. When g^p — > 0, the 
integrand is controlled by the perturbative functions in (9-10), which was also the case in 
§8. Contrastingly, when g NP — > oo, the integrand is controlled solely by the nonperturbative 
function, such that the integral of (9-10) is dominated by the b ?a region. Our task is 
to find a formula for finite g^p which interpolates between these two extreme cases, and, 
for this purpose, we perform the saddle-point evaluation of (9-10), assuming —gNpb 2 = 
—gNp(bl/Q 2 )e x ^ l3 ° as( ' Q2 ^ as a leading term in the exponent. As a result, we obtain (8-2), 
where C^(^) i s now replaced as 

C (o)( A ) - * MM + ?»r%e«rfm (9-11) 

and correspondingly the equation (8-5) to determine the saddle point is replaced by 

1-^ XsP = g^M e ^). (9-12) 

An important observation is that the ratio [h^(X) + Rn^X) + -Rat 2 (A)]/C ( - - ) (A) with 
(9-11) actually behaves as a quantity of order a s (Q 2 ) in the relevant region, A ps Xsp, of the 
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integration in (8-2), even for nonzero g^p — 0.5 GeV 2 . This fact supports the validity of 
the present approach, incorporating the nonperturbative function as (9-11). Also, this fact 
suggests that we may use an order counting that is similar to that explained below (8-7) 
for the contributions arising in the saddle-point evaluation. This is apparently correct when 
g^p is sufficiently small, but in general it has to be checked. Here we assume the same order 
counting as in §8, and check the accuracy of the resulting formula below, comparing it with 
the result of the numerical 6-integration of the original formula, (7-24) with (9-3) (see Table I 
below). 

Now we can perform the saddle-point evaluation of (8-2) at NLL accuracy and also a 
similar calculation starting from (7-38) for the unpolarized case. Immediately we find that 
the results (8-8)-(8T2) with (9-11) and (9-12) hold, up to the NNLL corrections corresponding 
to the 0(a s ) effects. 56 ** These results realize the "degree 0" approximation of §6.3. Also, it 
is straightforward to show that these results reproduce the results in the two extreme limits 
gNP and g^p — > oo, as desired for the "interpolation formula" mentioned below (9-10). 
Taking the ratio of the saddle point formulae (8-10) and (8-12), we obtain, 56 ) for the RHS of 
(9-9), 

^(0) \oz(26) ATjCNLh ^ Q2 ' y) 1 coz(26) SH ^^ %!M (9-13) 

TT " 2 COS(20) XNLL (0 , Q2, y) - 2 COS(20) H 6§/&| p ) ' (9 ^ 

up to NLL accuracy. Here, the common factor in (8-10) and (8-12), associated with the 
"very large perturbative effects" due to the Sudakov factor, as well as with the Gaussian 
smearing factor involving gNP, has canceled out. 

The new scale, bo/bsp, must be determined by solving (9-12) numerically, substituting 
A ( l ] = 2C F fr om (5-13) and the input values for Q and gNP, but it is useful to consider its 
general behavior: The LHS of (9-12) equals 1 at Xsp = 0, decreases as a concave function 
as Xsp increases, and vanishes at X$p = 1/[1 + / \2n (3 )] = 0.6, which corresponds to 
the solution (8-6) in the g^p = case, while the RHS is in general much smaller than 1 at 
\ SP = 0, increases as a convex function as Xsp increases, and is larger than 1 at Xsp — 1. 
Thus, the solution of (9-12) corresponds to the case with LHS = RHS ~ 1/2, more or less 



independently of the specific values of Q and gNP, so that we get bo/bsp — bo\/2g~NP- This 
result depends only weakly on the nonperturbative parameter g^p, and it suggests that one 
may always use b /bsp — 1 GeV for the cases of our interest, in which Q is on the order of 
several GeV and g NP ~ 0.5 GeV 2 , as in Figs. 2-5. The actual numerical solution of (9-12) 
verifies these simple considerations at the level of 20% accuracy. This fact is particularly 
helpful in the first attempt to compare (9-13) with the experimental data using (9-9), so as 
to extract the transversity distributions. 
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The above considerations also demonstrate that (9-10) is dominated by the region of b 
near b$p — 1 GeV ~ 1/tIqcd, corresponding to the boundary of perturbative and nonper- 
turbative physics. We can "safely" treat such a long-distance region, owing to the fact that 
the nonperturbative smearing (9-3) suppresses the too long-distance region, b ^> 1/Aqqb 
[see the discussion below (8-14)]. Moreover, the model dependence due to the specific choice 
of 9np largely cancels in the final result (9-13). 

Table I. The Q T -► limit of £?^ l {Qt) in (9-5) for RHIC and J-PARC kinematics. NB is 
obtained using the numerical 6-integration of (7-23), (7-24) and (7-38) with (9-3) and g^p = 0.5 
GeV 2 ; SP-I and SP-II are the results of the saddle-point formula (9T3) with (9T2) using the 
evolution operators from Q to bo/bsp, to NLL accuracy and to the customary NLO accuracy, 
respectively. 







Vs = 


200 GeV, y = 2 




Vs = 


: 10 GeV, 


y = o 


Q 


2GeV 


5GeV 


8GeV 


15GeV 


20GeV 


2GeV 


2.5GeV 


3.5GeV 


NB 


3.8% 


4.9% 


6.1% 


8.2% 


9.4% 


13.4% 


14.0% 


14.9% 


SP-I 


4.3% 


5.4% 


6.6% 


8.7% 


9.8% 


14.1% 


14.5% 


14.8% 


SP-II 


7.3% 


8.7% 


9.8% 


11.8% 


12.7% 


14.7% 


14.8% 


14.2% 



In Table I, "NB" lists ££jt L (<5t = 0) obtained from (9-5) using the numerical b- 
integration of (7-23), (7-24) and (7-38), with the kinematics of Fig. 4(a) for RHIC and of 
Fig. 5(b) for J-PARC. The result for Q = 5 GeV at RHIC and that for Q = 2 GeV at J-PARC 
coincide with the Qt — > limit of the dot-dashed curve in Figs. 3 and 5(a), respectively. 
The row labeled "SP-I" lists the results obtained using the saddle-point formula (9-13): We 
use bo/bsp obtained as the solution of (9-12) with g^p = 0.5 GeV 2 , and, according to (8-8), 
(8-9) and similar formulae for the unpolarized case, the parton distributions participating in 
(9T3) are obtained by evolving the customary NLO distributions at the scale Q to bo/bsp 
with the NLO evolution operators up to NLL accuracy, like (7-18). Here the "customary 
NLO distributions" are constructed as described above (9-3). Comparing NB and SP-I, we 
observe the remarkable accuracy of our simple analytic formula (9T3) for both RHIC and 
J-PARC, reproducing the results of NB to 10% accuracy, i.e., to the canonical size of 0(a s ) 
corrections associated with the NLL accuracy.*- 1 This supports the assumption concerning 
the order counting in our saddle-point evaluation, discussed below (9-12). 

Although the compact formula (9T3) is reminiscent of £^tt(Qt) i n (9- 7), which retains 
only the LL level resummation, or the QT-independent asymmetry of (9-8), it is different in 

*) If we use the fixed value bo/bsp = 1 GeV for all cases, instead of the solution of (9-12), the results in 
SP-I change by at most 5% for both RHIC and J-PARC kinematics. 
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the scale of the parton distributions from those leading-order results, and the only effect of 
this new scale bo/bsp allows SP-I to reproduce NB in Table I. Combined with this fact, the 
above derivation of (9-13), using (8-10) and (8-12), clearly demonstrates the characteristic 
features of the NLL soft gluon resummation effects on the asymmetries, which lead to the 
mechanism to "enhance" the dot-dashed curve in Figs. 3 and 5(a): In the Qt — > limit, 
the all-order soft-gluon-resummation effects on the asymmetry largely cancel between the 
numerator and the denominator of (9-13), but certain contributions at the NLL level survive 
the cancellation and are entirely absorbed into the unconventional scale bo/bsp — 1 GeV 
for the relevant distribution functions. Combined with the properties of ^tt L (Qt) in (9-9), 
this explains why ^tt L (Qt) is always larger than £&tt(Qt) i n (9 - 7) or Att in (9-8) as the 
NLL-level effects associated with the evolution [compare with the discussion below (9-7)]. 

Using (9-9), the saddle-point formula (9-13) can be directly compared with the experi- 
mental value of the asymmetries ^tt{Qt), observed around the peak of the Qt spectrum of 
the corresponding DY cross sections. This approach is particularly useful, because it does not 
require as large a numerical computations as that in §§9.1, 9.2. Furthermore, the approach 
with (9-13) reqires one point of caution about the accuracy, associated with the mismatch in 
classifying the terms between NLL and NLO. We recall that the parton distributions appear- 
ing in (9-13) are the NLO distributions up to the corrections at the NNLL level [see (8-8), 
(8-9)]. To obtain SP-I in Table I, the transversity distributions appearing in the numerator 
of (9-13) were obtained by evolving the customary NLO transversity Sq(x,Q 2 ) at the scale 
Q, to the scale b^/bsp using (7-18), which consists of the NLO evolution operators up to 
the NNLL corrections. Therefore, the formula (9-13) can be safely used in the region where 
the NNLL corrections are small; the NNLL corrections at Qt ~ correspond to 0(a s (Q 2 )) 
effects, applying the order counting of §§6.3 and 8, and should be negligible in general. How- 
ever, such a straightforward estimate might fail in the edge region of the phase space, as in 
the small x region. This follows from the fact that the relevant evolution operators (7-18) 
actually coincide with the leading contributions in the large-logarithmic expansion of the 
usual LO DGLAP evolution, *) (9-13) would not be accurate when the NLO corrections in 
the usual DGLAP evolution are large compared with the contributions of (7-18). Such a 
situation would typically occur in a region with small x\ 2 i corresponding to the case with 
large \/S. In Table I, we show other saddle-point results, "SP-II" , using the customary NLO 
distributions at the scale b^/bsp for the parton distributions participating in (9-13). Thus, 
SP-I and SP-II differ by contributions at the NNLL level. The results for SP-II indicate that 

*) This fact also suggests that one may use the fixed value bo/bsp — 1 GeV in (9-13) for all Q (and gNp) 
rather than solving (9-12) numerically for each different input value of Q and gNP, because the sensitivity 
of the LO evolution with respect to a small change of the scale is weak. 



72 



the NNLL corrections are moderate for large at RHIC, while they are expected to be 
small for small \fS at J-PARC. 

We propose that our simple formula (9-13) is applicable to analysis of low-energy exper- 
iments at J-PARC in order to extract the NLO transversity distributions directly from the 
data. On the other hand, (9-13) will not be very accurate for analyzing the data obtained 
from RHIC, but it will still be useful for obtaining the first estimate of the transversities. We 
emphasize that such a (moderate) uncertainty in applying our formula (9-13) to the RHIC 
case is not caused by the saddle-point evaluation, nor by considering the Qt — * limit, 
but rather is inherent in the general Qt resummation framework which, at the NLL level, 
implies the use of the evolution operators (7-18) with the LO anomalous dimension. A more 
accurate treatment of the small-x region of the parton distributions relevant to the RHIC 
case would require the resummation formula to NNLL accuracy, where the NLO anomalous 
dimensions participate in the evolution operators (7-18) from Q to b /b (see e.g. Ref. 35)). 



§10. Conclusions 



The era of "testing QCD" is clearly finished and today the improvement of the precision 
of predictions is one of the most important subjects in QCD. We have reported our analysis 
of the transversely polarized DY process, which can be measured in ongoing experiments at 
RHIC and in future experiments at J-PARC, GSI, etc. We have demonstrated that the QCD 
corrections relevant for the tDY at a measured Qt are now under control over the entire 
range of Qt, so that the dilepton Qt spectrum as well as the Qt dependence of the double- 
spin asymmetry for the tDY is predictable in a systematic QCD factorization framework. In 
particular, in the small Qt region, where most of the dileptons are produced in experiments, 
these observables are dominated by universal logarithmically-enhanced contributions, which 
we have taken into account by an all-order resummation of multiple soft-gluon emission. 
Based on this framework, we have performed a detailed numerical study of the observables 
in the tDY at RHIC and J-PARC, using a model of the transversity which is constructed 
to satisfy a known theoretical constraint. We find that the logarithmically-enhanced higher- 
order QCD corrections associated with soft gluon resummation drive a mechanism that 
amplifies the double transverse-spin asymmetries in the small Qt region. Quantitatively, 
our results indicate that the double transverse-spin asymmetries seem to be small at RHIC, 
while they will be large at J-PARC. This would imply that the measurement of the tDY 
double transverse-spin asymmetries at RHIC might require considerable effort. Still, we 
expect the data from RHIC, as well as that from J-PARC, GSI, etc., that are comparable 
with our predictions. In particular, we have proposed a remarkably simple formula for 
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the double transverse-spin asymmetries, which provides us with a new direct approach to 
extracting the transversity distributions from experimental data. We believe that the results 
presented in this paper will be useful not only for experiments at RHIC and J-PARC but 
also for future precision measurements. 

We have obtained all our results concerning the QCD corrections for the DY in the MS 
scheme, for the first time, directly in the dimensional regularization with transverse polariza- 
tion, and our calculation technique can be straightforwardly applied to other polarized and 
unpolarized processes. As described in this paper, the formalism of transverse-momentum 
resummation we employed is the same for many other polarized and unpolarized processes, 
and has a wide range of application, in particular, to the processes observed at LHC, where 
rigorous treatments of multi-scale problems in QCD are mandatory and require resumming 
ubiquitous soft gluons. QCD is now facing a new challenging era regarding its precision and 
applicability. We hope that various kinds of new experiments and theoretical investigations 
will be able to clarify not only perturbative and nonperturbative aspects of QCD but also 
the full structure of all interactions in Nature. 
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Appendix A 

The Transversity Distribution 

The transversity distribution of a quark inside a spin- 1/2 hadron h is defined in terms of 
a matrix element of a non-local light-cone operator as 14 ) 

5q h (x,f, 2 ) = J ^e iX %h{P,S)\m[^Mil,m^n) \h(P,S)) , (A-l) 

where P and S are the momentum and transverse spin of the hadron h, respectively, and n M 
is a light-like vector which satisfies n 2 = 1 and P ■ n = 1; for example, when we ignore the 
nucleon mass, setting P 2 = 0, and work in a frame where P» = (P+/ v / 2)(l, 0, 1), we have 
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S 11 = (0, 5,0) and = (l/y/2P + )(l,0, -1). The path-ordered link operator (Wilson line) 

[0, \n] = V exp (-ig J d\'n ■ A(\'n)^j , (A-2) 

which connects the two quark operators -0(0) and ip(\n), is required to ensure the gauge 
invariance of the non-local operator. In the definition (A-l), the scale fu, 2 of the distribution 
corresponds to the renormalization scale of the non-local operator on the RHS. 

The "quark transversity-distribution inside a quark" can be obtained by evaluating the 
matrix element (A-l), with h = quark q', in perturbation theory with an appropriate reg- 
ularization of the IR divergence. Employing the dimensional regularization in D = 4 — 2e 
dimensions for the IR and UV divergences and evaluating the matrix element up to one-loop 
accuracy, (A-l) for h — q' yields*) 

5q q >(x,[i 2 ) = 5 qq > 

in the MS scheme for renormalization of the non-local operator of (A-l), where 1/e = 1/e — 
7s + ln(47r), and A T P qq (x) is the LO splitting function for the transversity: 

A T P qq {x) = C F [j^U + I 5(1 ~ X) ] ■ (A ' 4) 

We note that we can derive the DGLAP evolution equation at the LO level by differentiating 
(A-3) with respect to ji 2 . The result is [compare with (5-15)] 

V 2 ^5 qh (x^ 2 ) = ^1 f X ^A T P qq (f) 6q h (z,^). (A-5) 

Together with the density distributions (qh(x,fj, 2 ), gh(x,fi 2 )) and the helicity distribu- 
tions (Aqh(x, fi 2 ), Agh(x, fi 2 )), the transversity distribution 5qh(x,fi 2 ) forms a complete 
set of twist-2 parton distributions for a spin- 1/2 hadron h. Here the quark distributions 
qh(x,fi 2 ) and Aqh(x,fi 2 ) are defined by the matrix element (A-l) for the unpolarized and 
longitudinally-polarized hadron states, with the chiral-odd Dirac structure ^75$ replaced 
by the corresponding chiral-even structures i/i and ^75, respectively. Thus, unlike these 
two chiral-even quark distributions, the transversity is chiral-odd. Because all gluon distri- 
butions at twist-2 are chiral-even, and the QCD interaction conserves chirality, the gluon 
distributions do not mix into the DGLAP evolution (A-5) for the transversity, even at NLO 
or higher-order level. As a result, the QCD evolution of the transversity beyond the LO is 
of the non-singlet type; i.e., it involves only the mixing between the quark and antiquark 

*) A similar calculation for the unpolarized quark distribution q q /(x,[i 2 ) is described in Ref. 1). 
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distributions associated with the same flavor, and it is much simpler than the evolution for 
the other two quark distributions. For example, in the Mellin iV-space, the NLO evolution 
of the transversity distribution is written as 
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where Sq^ifJ 2 ) = dxx N 1 5q ± (x, fx 2 ) with 5q ± (x,fj, 2 ) = 5qh(x, /i 2 ) ± 5qh(x, fu, 2 
anomalous dimension is given by 70 ) 
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where r\ = ±, and the harmonic sums are defined by 
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The first moment of the LO splitting function is negative (^ q °J N=1 = — |), so that the 
transverse polarization in total (the "tensor charge") decreases as the scale increases.*-* 



(A-6) 
The LO 
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*) A pedagogical explanation of the transverse-spin flip by collinear radiation is given in the appendix 
of Ref. 72). 
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We note that the dominant large- z behavior (5-17) of the DGLAP splitting functions 
(5-16) reflects the dominant large N behavior of the corresponding anomalous dimensions 
(7-13). For the present case, the anomalous dimensions given in (A-7) and (A-8) at large N 
behave as 



? (D 

lqq,N - K 1111 " T W) ^ 

where the ellipses stand for the terms that vanish as N — > oo, and 



A = Q*N + lE ) - -J- + • • • , (A- 12) 



7^ = -4 2) ( lnAr + 7£) + ---, (A-13) 

where the ellipses stand for the corrections of 0(1). The coefficients on the RHS of (A-12) and 
(A-13) coincide with (5-13). The same asymptotic behavior arises in the anomalous dimen- 
sions for the density and helicity distributions q h (x, p 2 ) and Aq h (x, p 2 ) (see e.g. Ref. 73)). 

Appendix B 

The Cross Section Integrated over the Rapidity y 



Integrating (4-29) over the rapidity y of the dileption in the final state, we obtain the 
rapidity-integrated cross section in the MS factorization scheme, 

dQHQ%d<p k =^ cos ( 2 <fc-^-^) Jdy[A T X(Q 2 T ,Q 2 ,y) + A T Y(Q 2 T ,Q 2 ,y)}, (B-l) 
including the 0{a s ) QCD mechanism. We use the following integration formulae: 

dy5H{x\ ) x\) = J dx\dx 2 SH(xi, x 2 ) 8{x\x 2 — r), (B-2) 
/ dy f -A T P qq (z)SH(^,x{] = [ ^lA T P qq (r 12 )6H(x 1 ,x 2 )9(x 1 x 2 -r), 

J J x o Z \ Z J J X\X 2 

(B-3) 

f dy f 1 — T ^H{x l ,x* 2 ) _ I dx x dx 2 t 12 SH(xi,x 2 ) q(^. x _ T ) 
J J^fTev xi - xf xix* 2 J xix 2 a/(1 - r 12 ) 2 - 4pi2 1 2 

(B-4) 
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where 5H(xi,x 2 ) denotes the product of the transversity distributions, (2-4), suppressing 
the scale dependence, 5Hi is defined in (4-25), r is given in (4-7), and we have 
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We thus get 
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For Q T > 0, the delta function 5{Q T ) involved in (B-7) vanishes, and, after some calculation, 
(B-l) with (B-7) and (B-8) reduces to 
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This coincides with the result obtained by imposing the restriction Q\ > in Ref. 16), and 
it provides the LO QCD prediction in the large Qt region. 

To obtain the QCD prediction for the rapidity-integrated tDY cross section over the entire 
region of Qt-, we integrate (7-25) over y. This gives the rapidity-integrated NLL+LO cross 
section in the MS scheme. The y- integral of A T X NLL (Q T , Q 2 ,y), as well as A T Y(Q T , Q 2 ,y), 
can be performed straightforwardly using (B-2)-(B-5), similarly to the above fixed-order 
case. For AtX nll (Q T ,Q 2 ,y), the result can be expressed in the Mellin moment space 
with respect to t — x\x\ at fixed Q: Changing the integration variables from x\ and x® 
to t and y in the definition of the double Mellin moment (7-23) [see (7-6)], we obtain 
/ dx\dxl(x\) N ^\xl) N *- l A T X^ h (Q 2 T , Q 2 , y) = f drdy t (n 1+ n 2 )/2-i e (N^N 2 ) y ^nll^ ? g 2? y y 
Setting Ni = N 2 in this relation, we get the Mellin-space representation of the corresponding 
y-integral as 

^dTT N ^\J dyA T X^ L (Q 2 T ,Q 2 ,y) 
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Jo 



1 + ^^lc F (ir 2 - 8) 



6H NuNl (Q 2 )I NuNl (Q 2 T ,Q 2 ), (B-10) 



withI NuNl (Q 2 T ,Q 2 ) of (7-24). 



Appendix C 

Cancellation at Q T = in Eq. (7-26) 

In this appendix, we show that the terms proportional to the delta function S(Q^) cancel 
out in the difference A T X — Z\ T X NLL | fo in (7-26), and as a result, A T Y differs from A T Y 
only by the terms that are less singular than 1/Q% or 5(Q\) as Q\ — > 0. This shows that the 
NLL+LO cross section (7-25) is smooth and well-defined in the entire region of Qt-, including 
Qt = 0. For this purpose, we explicitly show that [see (7-27)] the relation 

dQ 2 T (I n - X n ) = lim f - J 2 dQ T T n + J 2 dQ T T n J = (C-l) 

holds for n — 1 and 2, which means that the delta function 5(Q T ) implicit in the + distri- 
butions in X12 cancels with that in Xi j2 . 

The first equality in (C-l) follows from (7-33). To show the second equality, it is conve- 
nient to consider the generating functions of 1 n and X n . The explicit form of the former is 
already given in (6-4) as 
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and r(l + ?7)/r(l - 77) = (6§/4)" exp [-2 ^g^V r+1 ] • The Q T integration of (C-2) 
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and the Qr integration of this function is obtained as 
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Taylor expanding (C-4) and (C-8) about 77 = and comparing the results, we see that (C-l) 
is valid for n — 1 and 2. For n > 3, (C-l) is not valid, because terms including the zeta 
function ((n) appear differently in (C-4) and (C-8). 
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